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A theory  and  computational  treatment  of  photodesorption  of  neutral  molecules  from 
metal  surfaces  caused  by  absorption  of  femtosecond  visible  or  UV  light  pulses  are 
presented  in  this  work.  The  molecule  and  metal  system  are  separated  into  primary  and 
secondary  regions.  This  theory  introduces  the  Liouville-von  Neuman  equation  for  the 
density  operator  and  derives  from  it  the  effective  equations  for  the  two  coupled  regions, 
assuming  a self-consistent  factorization  of  the  density  operator.  The  secondary  region 
is  described  as  a stochastic  medium  providing  energy  dissipation  during  the  desorption. 
Direct  and  indirect  desorption  mechanisms  are  considered,  with  the  latter  derived  from 
a local  electromagnetic  field  which  excites  the  metal  substrate  eleetron-hole  excitations. 
These  couple  through  Coulomb  interactions  to  the  primary  region,  to  produce  an  effeetive 
transition  dipole,  following  which  the  adsorbate  bond  is  broken.  We  derived  equations 
for  this  effective  primary  dipole  by  applying  first  order  perturbation  theory  for  small 
fields  and  by  including  all  orders  for  larger  fields. 

The  primary  time-dependent  Schrodinger  equation  was  numerieally  solved  by  prop- 
agation of  wavepackets,  using  the  split-operator  propagation  scheme  to  construct  the 
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time  evolution  operator,  with  a modification  to  include  imaginary  dissipative  potentials 
and  a time  dependent  interaction  with  light.  The  diabatic  electronic  representation  was 
used  in  the  time  propagation.  A Fortran  77  code  for  the  wavepacket  propagation  in 
two  degrees  of  freedom  over  two  coupled  complex  PESs  was  developed  and  applied 
to  the  case  of  CO  desorbing  from  the  Cu(OOl)  surface.  This  code  allows  for  a quan- 
tum mechanical  calculation  of  any  two  dimensional  desorption  process  with  Cartesian 
coordinates.  A numerical  algorithm  which  combines  the  split  operator  scheme  with 
a discrete  variable  representation  is  also  developed  for  the  rotational  motion  of  the 
desorbing  species. 

As  applications,  the  calculations  of  the  CO  photodesorption  were  done  varying 
several  different  parameters:  the  light  pulse  fluence,  width,  and  central  frequency,  the 
magnitude  of  the  dissipative  imaginary  potentials  and  their  coupling,  and  the  temperature 
of  the  initial  system  and  the  isotope  masses  of  CO.  Results  of  the  calculations  are  in 
good  agreement  with  experimental  data  for  photodesorption  yield  dependence  on  light 
pulse  fluence,  femtosecond  response  time  and  for  the  anomalous  photodesorption  yield 
dependence  on  CO  isotope  masses.  The  photodesorption  yield  dependence  on  light  pulse 
chirping  is  also  given  and  shows  that  positive  and  negative  chirping  have  pronounced 
effects  on  CO  yields,  with  positive  chirping  decreasing  and  negative  chirping  increasing 
the  yields,  respectively. 
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CHAPTER  1 
INTRODUCTION 

Photostimulated  desorption  processes  of  small  molecules  from  well-defined  sub- 
strate metal  surfaces  have  attracted  intense  interest  in  recent  years  due  to  interest  on 
elementary  chemical  processes  at  surfaces  on  the  ultrafast  timescale  of  femtoseconds, 
and  to  provide  complementary  insights  into  molecule-surface  interaction  dynamics.  De- 
spite considerable  experimental  and  theoretical  effort  [1-7],  a complete  understanding 
of  these  processes  remains  elusive.  Recent  femtosecond  photodesorption  experiments 
have  shown  some  general  new  properties  including  a nonlinear  dependence  of  desorp- 
tion yield  on  laser  fluence,  significantly  enhanced  desorption  yield  over  that  obtained 
with  nanosecond  laser  pulses  with  comparable  fluence,  and  a femtosecond  desorption 
response  time.  These  observations  are  inconsistent  with  a conventional  thermal  or  pho- 
tochemical mechanism  for  desorption  with  picosecond  or  nanosecond  laser  pulses.  The 
purpose  of  this  study  is  to  understand  the  femtosecond  light  pulse  photodesorption  of  a 
molecule  bonded  to  a metal  surface  and  to  provide  computational  methods  for  related 
physical  properties. 

System  and  Phenomena  of  interest 

In  this  dissertation  we  investigate  the  desorption  induced  by  visible  or  UV  fem- 
tosecond light  pulses,  in  which  case  the  system  is  electronically  excited.  Indicating  with 
M the  absorbed  molecule  and  with  S the  metal  substrate,  the  relevant  processes  are 


1 


2 


M/5'  + ^ > M* / S {photon  absorption) 

M* I S S + M {desorption) 


^M^jS  {el  — ^ vih  transfer) 
M^S-^M/S  { vihra.tiona.1  rela.xo.ti.on ) 


—^S  + M {desorption) 

where  (/>  is  a visible  or  UV  photon,  and  we  show  the  possibility  of  electronic  to 
vibrational  energy  transfer  before  desorption. 

The  process  of  interest  is  the  indirect  mechanism  which  has  been  observed  experi- 
mentally for  photons  off  resonance  with  the  bond  excitation  energies  in  CO/Cu(001).  In 
our  model,  the  CO  molecule  is  absorbed  at  an  on-top  location  perpendicular  to  the  sur- 
face. The  electronic  states  of  the  whole  system  can  be  understood  by  coupling  the  a and 
7T  molecular  orbitals  of  CO  and  the  3d/4s  and  4p  electronic  energy  bands  of  the  metal. 
In  principle  there  is  also  a coupling  between  vibrations  of  CO  and  of  the  metal  lattice, 
but  we  shall  argue,  as  suggested  by  experiments,  that  metal  vibrations  are  too  slow  to 
participate  in  the  desorption  process,  and  we  shall  ignore  them  in  the  present  model. 


Overview  of  Previous  Work 

The  theory  of  photodesorption  has  been  developed  by  several  authors  over  the 
years,  in  connection  with  experiments  on  desorption  induced  by  electronic  transitions 
(DIET)  and  desorption  induced  by  multiple  electronic  transitions  (DIMET)[l-7].  Basic 
questions  arising  in  these  studies  [8—27]  have  been  whether  the  excitations  of  the 
substrate  are  predominantly  vibrational  or  electronic;  whether  desorption  occurs  directly 
at  the  adsorbate  site  or  indirectly  with  electronic  excitation  occurring  first  in  the  substrate; 
and  what  is  the  nature  of  the  electronically  excited  states  leading  to  desorption.  These 
aspects  are  being  probed  by  experiments  with  light  pulses,  where  time  scales  become 
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apparent  and  the  vibrational  and  electronic  excitation  mechanisms  can  be  separated.  It 
has  recently  become  clear  that  for  bond  breaking,  a fast  event,  electronic  excitations 
are  the  dominant  ones  [28-39]. 

Unlike  conventional  photochemistry  using  longer  pulses  and  lower  intensities 
(DIET)  which  lead  to  linear  dependence  of  the  yield  and  independence  of  the  final 
state  distributions  on  laser  fluence  [8-11],  the  desorption  induced  by  femtosecond  laser 
pulses  (DIMET)  has  been  studied  for  several  systems  [28-32,  34,  37-39],  and  has  led  to 
the  demonstration  of  novel  desorption  dynamics  for  molecules  on  metal  surfaces.  Some 
properties  have  been  observed  in  these  studies:  a nonlinear  dependence  of  desorption 
yield  (T)  on  laser  fluence  (F),  Y (x  F''  with  v reported  in  the  range  3 to  8.4,  and  a 
very  short  system  response  time,  300  fs  to  1.8  ps;  and  significantly  enhanced  desorption 
yield  over  that  obtained  with  nanosecond  laser  pulses  with  comparable  fluence. 

The  first  direct  femtosecond  measurement  of  the  time  evolution  of  a surface  chemical 
reaction,  the  making  and  breaking  of  chemical  bonds  on  metal  surface  was  reported  by 
Prybyla  et  al.  [29].  They  used  second-harmonic  generation  as  a time-resolved  probe 
of  the  coverage  of  absorbed  molecules.  They  were  thus  able  to  measure  the  time 
of  desorption  with  a resolution  limited  only  by  the  accompanying  second-harmonic 
response  due  to  substrate  heating.  The  experiments  were  performed  on  the  CO/Cu(l  1 1) 
system  at  coverage  of  0.52  monolayer(ML),  at  a temperature  of  95  K,  with  100  fs 
laser  pulses  operating  at  around  620  nm.  They  found  that  CO  molecules  desorb  from  a 
Cu(lll)  surface  within  325  fs  after  irradiation  with  those  100  fs,  2 eV  photon  energy 
pulses.  They  also  found  that  the  desorption  yield  measured  with  a mass  spectrometer 
increase  nonlinearly  with  the  intensity  and  could  be  fitted  by  the  power  law  I For 
the  case  of  a fluence  of  4.5  mJ/cm^,  the  yield  is  0.05  ML. 
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Other  femtosecond  time-resolved  measurements  of  desorption  performed  by  differ- 
ent groups  on  different  system  include  NO/Pd(lll)  [28,  30,  37],  02/Pt(l  1 1)[34,  38], 
02/Pd(lll)  [32]  and  CO/Cu(100)  [39].  All  these  experiments  demonstrate  similar  des- 
orption dynamics  as  given  in  the  work  by  Prybyla  et  al.[29]. 

Several  models  have  been  introduced  to  describe  photodesorption  induced  by  light 
pulses.  A model  based  on  excited  hot  electrons  scattered  into  an  unoccupied  orbital 
of  the  adsorbate  assumes  the  formation  of  a temporary  negative  ion  and  develops  a 
semiclassical  wavepacket  dynamics  for  the  subsequent  desorption  [12-14].  A model 
assuming  multiple  electronic  transitions  is  based  on  the  coupling  of  the  center  of  mass 
of  the  adsorbate  with  the  electron-hole  excitations  of  the  substrate  and  introduces  a 
friction  derived  from  electron  phase  shifts  at  the  Fermi  energy  of  the  metal;  this  is  used 
in  a Langevin  equation  which  is  solved  numerically  [15-17].  Femtosecond  desorption 
has  been  described  with  trajectory  calculations  incorporating  electronic  friction  and 
fluctuation  forces,  where  the  friction  is  calculated  with  a molecular  orbital  procedure 
[20-22].  An  alternative  approach  using  the  density  operator,  which  has  grown  from  the 
interpretation  of  experiments  in  terms  of  substrate  temperatures  for  electrons  and  the 
lattice,  also  introduces  dissipation  in  a similar  fashion  as  we  will  explain,  but  with  a 
different  treatment  where  a phenomenological  time  dependent  dissipative  term  is  added 
in  the  density  matrix  equations  [24,  25]  instead  of  deriving  it  directly  from  the  system 
Hamiltonian. 

In  what  follows  a theory  is  developed  for  the  dynamics  of  photodesorption  by  light 
pulses,  showing  how  relevant  properties  can  be  calculated,  and  present  a model  where  S 
is  the  surface  Cu(OOl),  and  M is  the  CO  diatomic.  The  approach  incorporates  some  of 
the  same  aspects  introduced  in  studies  of  the  photodissociation  of  polyatomic  molecules 
[40^2].  The  system  is  separated  into  a primary  region  (the  adsorbate  and  bonding 
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substrate)  and  a secondary  region  (the  remaining  substrate)  which  are  strongly  coupled 
and  interacting  in  a selfconsistent  way.  We  concentrate  on  the  dynamics  of  the  primary 
region,  governed  by  an  effective  Hamiltonian  which  includes  dissipation  and  can  be 
derived  under  the  assumption  that  the  secondary  region  undergoes  stochastic  motions 
over  which  it  is  possible  to  average  primary  quantities  [43]. 

Properties  of  interest  are  the  photodesorption  yields  by  light  pulses,  the  survival  and 
transition  probabilities,  electronic  distributions  on  electron  states,  and  the  distribution 
of  kinetic  and  vibrational  energies  of  the  desorbed  molecule.  Some  of  these  properties 
have  been  calculated  by  us  for  an  assumed  direct  mechanism  of  desorption  of  CO 
from  the  Ni(OOl)  surface  by  photons  with  energies  around  the  C — Ni  bond  excitation 
energies[44],  and  indirect  mechanism  of  CO  from  the  Cu(OOl)  surface  by  photons  with 
energies  off  C — Cu  bond  excitations  energies  [45]. 

In  this  study,  a new  approach  is  developed  for  indirect  desorption  by  introducing 
an  effective  primary  transition  dipole,  which  comes  from  the  coupling  of  the  primary  or 
secondary  regions.  A detailed  theoretical  treatment  of  the  secondary  region  electron-hole 
excitation  on  the  effective  primary  transition  dipole  is  given,  by  treating  the  interaction 
between  the  light  pulse  and  the  secondary  region  with  both  first  order  perturbation  the- 
ory and  nonlinear  response  theory.  The  effective  primary  region  Hamiltonian  including 
dissipation  in  the  secondary  region  and  the  indirect  light  pulse  effect  on  the  primary 
region  is  derived  and  applied  to  simulate  the  CO/Cu(001)  femtosecond  photodesorption 
process.  The  details  of  the  electronic  structure  of  both  ground  and  excited  electronic 
states  of  CO/Cu(001)  have  been  extracted  from  calculations  using  semi-empirical  molec- 
ular orbital  theory  [46-48]  in  which  the  metal  substrate  is  described  by  a cluster  of  Cu 
metal  atoms.  Calculations  of  the  quantum  dynamics  of  desorption  have  been  carried 
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out  by  propagating  wavepackets  using  the  split-operator  scheme  [49,  45]  with  some 
modifications. 


Outline  of  the  Thesis 

Chapter  2 describes  the  general  physical  model  for  adsorbate  molecule  photodes- 
orption from  a metal  surface  or  an  extended  system.  The  general  treatment  is  based 
on  separating  the  extended  system  into  primary  region  (the  adsorbate  and  bonding  sub- 
strate) and  secondary  region  (the  remaining  substrate).  By  introducing  the  Liouville-von 
Neumann  equation  for  the  density  operator  and  assuming  a self-consistent  factorization 
of  the  density  operator,  the  effective  equations  for  the  strongly  coupled  primary  and  sec- 
ondary regions  are  derived.  The  secondary  region  is  described  by  a stochastic  medium 
which  provides  energy  coupling  and  dissipation  during  the  desorption.  Two  mechanisms 
of  photodesorption  are  discussed,  the  direct  and  indirect  mechanism.  In  the  former,  the 
photons  are  assumed  to  be  absorbed  directly  into  the  primary  region  and  excite  it,  while 
in  the  latter,  the  photons  are  absorbed  by  the  secondary  region,  and  excite  the  substrate 
electron-hole  excitations,  following  which  the  adsorbate  bond  is  broken. 

Chapter  3 presents  the  dynamics  of  photodesorption  for  the  primary  region.  The 
primary  region  time-dependent  Schrodinger  equation  with  an  effective  time  dependent 
Hamiltonian,  including  dissipation  from  the  secondary  region  and  interaction  with  a 
light  field,  is  described  within  electronic  adiabatic  and  diabatic  representations  for  a 
system  with  two  potential  energy  surfaces.  The  expressions  for  survival  probabilities, 
population  on  both  potential  energy  surfaces  and  photodesorption  yields  are  presented. 

Chapter  4 presents  the  general  theory  of  the  secondary  region  interaction  with  both 
the  primary  region  and  the  time-dependent  light  field.  Describing  the  electron-hole 
excitations  by  a time-dependent  one-electron  density  operator,  the  dissipative  term  of 
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the  primary  region  Hamiltonian  has  been  formally  derived  by  using  stochastic  dynamics 
treatment.  The  effective  transition  dipole  for  the  primary  region  is  also  introduced  by 
developing  both  linear  response  theory  and  nonlinear  response  theory  of  the  secondary 
region  to  the  light  field.  In  both  cases,  the  effective  transition  dipole  was  found  to  be 
either  enhanced  from  the  original  transition  dipole  or  became  nolinearly  dependent  on 
light  pulse  fluence. 

Chapter  5 presents  the  computational  aspect  of  this  work.  A modified  split- 
operator  scheme  is  developed  for  the  numerical  time  evolution  of  wavepackets  with  time- 
dependent  Hamiltonian  in  Cartesian  coordinates.  For  problems  with  angular  degrees  of 
freedom,  an  algorithm  has  been  introduced  which  combines  the  split-operator  method 
with  a discrete  variable  representation.  A short  description  is  also  given  on  how  to 
implement  the  modified  split-operator  scheme,  and  the  flow  chart  of  a Fortran  77 
program  is  given  in  Appendix.  In  this  chapter,  I also  give  the  initial  conditions  for 
the  propagation  with  the  physical  implications. 

Chapter  6 presents  an  application  of  the  method  to  the  indirect  photodesorption  of 
CO  from  Cu(OOl).  A semiempirical  parametrization  of  the  complex  potential  energy 
surfaces  is  presented.  The  real  parts  were  obtained  from  the  ZINDO  electronic  structure 
calculations  for  a COCu  molecule  and  CO  absorbed  on  a Cug  cluster  [45].  The  results 
were  then  rescaled  to  fit  the  known  experimental  results.  For  the  imaginary  parts  and  the 
effective  transition  dipole,  theoretical  models  were  introduced  based  on  previous  chapters 
in  this  study  and  were  parameterized  with  known  experimental  results.  Numerical  results 
were  then  presented  and  included  properties  such  as  photodesoption  yields,  population  on 
either  electronic  state  and  time  evolution  of  the  amplitude.  Results  are  also  presented  for 
properties  dependent  on  the  pulse  fluence,  width,  frequency,  chirping  and  dissipation 
strength. 
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Chapter  7 presents  the  results  on  the  dependence  on  CO  isotope  masses  and  initial 
substrate  temperature.  The  results  show  good  agreement  with  preliminary  experimental 
results  which  show  anomalous  desorption  yield  dependences  on  CO  isotope  masses. 

Chapter  8 discusses  advantages  and  distinctive  features  of  the  present  approach,  and 
describes  possible  improvements  of  this  work. 
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CHAPTER  2 

PHYSICAL  MODEL  FOR  PHOTODESORPTION 

This  chapter  overviews  the  general  treatment  for  adsorbate  molecule  photodesorption 
from  a metal  surface  or  an  extended  system[40-43,  45]  and  generalizes  it  to  explicitly 
include  light  field  driven  terms  which  are  necessary  to  treat  femtosecond  photodesorption 
dynamics. 

The  general  treatment  is  based  on  separating  the  extended  system  into  primary 
and  secondary  regions.  By  introducing  the  Liouville-von  Neumann  equation  for  the 
density  operator  and  assuming  a self-consistent  factorization  of  the  density  operator, 
the  effective  equations  for  the  strongly  coupled  primary  and  secondary  regions  are 
derived.  The  secondary  region  is  described  by  a stochastic  medium  which  provides 
energy  coupling  and  dissipation  during  the  desorption.  There  are  two  mechanisms  by 
which  photodesorption  can  occur,  the  direct  and  indirect  mechanism.  In  the  former, 
the  photons  are  assumed  to  be  absorbed  direetly  by  primary  region  and  excite  primary 
region,  while  in  the  latter,  the  photons  are  absorbed  by  the  secondary  region  and  excite 
the  substrate  electron-hole  excitations,  following  which  the  adsorbate  bond  is  broken. 

Partition  of  the  System  into  Primary  and  Secondary  Regions 

The  basic  model  in  this  study  consists  of  an  adsorbate,  a substrate  and  a light 
field.  Depending  on  whether  the  mechanism  is  direct  or  indirect,  either  the  adsorbate 
is  simultaneously  coupled  to  the  substrate  and  the  light  field,  or  the  substrate  is 
simultaneously  coupled  to  the  adsorbate  and  the  light  field.  In  both  case,  since  the 
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light  field  is  strong,  it  is  sufficient  to  treat  the  light  field  semiclassically.  The  models  are 
schematically  illustrated  in  Fig.  (1)  and  Fig.  (2)  for  direct  and  indirect  photodesorption 
respectively. 

The  total  Hamiltonian  for  the  extended  system  including  the  light  interactions  with 
both  the  adsorbate  and  the  substrate  can  be  written  in  five  terms. 


/N 


/s  /N  y\ 

M.  r -L  I 


/N  /N  /N  /S  /\ 

where  Hm,  Hs,  H\fs{t),  H\n{t)  and  Hsiit)  are  the  Hamiltonians  for  the  free 
adsorbate  molecule,  the  free  substrate  surface,  the  coupling  between  the  adsorbate  and 
substrate  surface,  the  coupling  between  the  adsorbate  and  the  light  field,  and  the  coupling 
between  the  substrate  surface  and  the  light  field.  Instead  of  working  with  this  expression, 
it  is  more  convenient  to  rearrange  the  terms  to  reflect  the  fact  that  the  absorbed  molecule 
is  strongly  coupled  (forming  a chemical  bond)  to  the  adjacent  substrate.  This  region 
changes  extensively  during  desorption  and  requires  a detailed  treatment  in  terms  of 
local  degrees  of  freedom,  while  the  more  distant  regions  of  the  metal  undergo  collective 
excitations.  Introducing  a primary  (p)  region  which  includes  the  adsorbate  and  the 
adjacent  substrate,  and  a secondary  (s)  region  including  the  remaining  substrate,  the 
total  Hamiltonian  become 


H[t)  — Hp  + + Hpg[t)  + Hp{t)  + Hg[t)  (3) 

where  now  the  first  term  is  the  Hamiltonian  of  the  p-region,  the  second  term  is  the 
Hamiltonian  for  the  s-region,  the  third  term  is  the  coupling  between  the  primary  and 
secondary  regions,  the  fourth  term  is  the  coupling  between  the  primary  region  and  the 
light  field,  and  the  last  term  is  the  coupling  between  the  secondary  region  and  the  light 


field. 
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The  state  of  the  extended  system,  consisting  of  the  p-  and  s-  regions  can  be  described 

/\ 

by  the  density  operator  T.  The  dynamics  of  T is  governed  by  the  Liouville-von  Neumann 
(L-vN)  equation 


= f(^)] 


(where  h = 1 here  and  henceforth)  which  must  be  solved  with  initial  conditions 
corresponding  to  thermal  equilibrium  at  time  t = 0 before  the  arrival  of  a light  pulse, 
that  is 


r(0)  - Te,  (5) 

The  coupling  of  the  p-  and  s-regions  can  not  be  assumed  to  be  weak.  Each  region 
drives  the  other.  Energy  dissipation  in  the  s-region  can  only  be  properly  described  by 
recognizing  that  the  p-region  excitation  forces  couplings  among  the  s-region  collective 
motions,  which  in  turn  leads  to  energy  dissipation  in  that  region;  the  energy  transfer  to 
the  p-region  due  to  the  s-region  excitation  by  the  external  light  field  must  be  described 
by  coupling  the  s-region  electron-hole  excitation  with  the  motion  of  the  p-region,  and 
this  leads  to  energy  transfer  to  that  region.  A non-perturbative  treatment  can  be  based 
on  a molecular  self-consistent  field  (SCF)  procedure  which  factors  the  density  operator 
into  p-  and  s-  parts, 

f(f)  = P{t)  (g)  (6) 

with  each  factor  satisfying  an  L-vN  equation  with  an  effective  Hamiltonian  obtained 
from  a variational  treatment  [41,  43].  The  phase  factor  including  x can  be  expressed 
as  an  integral  over  time  but  will  cancel  in  our  expressions  for  properties,  and  need  not 
be  considered.  For  the  p-region  this  gives 
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and  for  the  s-region 


/s 


The  effective  Hamiltonian  for  the  p-region  is  given  by. 


+ Gp{t) 


where  Gp{t)  comes  from  the  coupling  between  the  s-region  and  the  p-region.  It  gives 
the  energy  coupling  and  dissipation  effects  of  the  s-region  on  p-region. 

For  the  density  operator  we  separate  the  effects  of  the  external  light  field 

from  the  p-  and  s-  region  coupling,  as 


r'(()  = rs(()  + r,"(() 


where  rQ(f)  represents  the  s-region  density  operator  in  the  absence  of  the  light  field,  and 

yN 

r*(/)  represents  the  s-region  density  operator  changes  due  to  the  light  field  interaction 
with  the  extended  system  causing  changes  in  the  electronic  density  distribution.  The 
effective  Hamiltonian  for  the  p-region  can  now  be  written  as, 

Fp{t)  = Hp  + Hp{t)  -h  Gp{t)  -f  Gp{t) 

Gj(()  = (r.[f5(()//„,(«)|/«r,[f '(()]  (1 1) 

where  Gp{t)  is  due  to  the  coupling  between  the  s-region  and  the  p-region  in  the  absence 
of  the  light  field  and  G^,{t)  arises  from  the  effects  on  the  p-region  of  the  interaction  of 
the  light  field  with  the  extended  system. 
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Similarly,  the  effective  Hamiltonian  for  the  s-region  is  given  by. 


F,{t)  = H.,  + Hl{t)  + G,{t) 


where  comes  from  the  coupling  between  the  s-region  and  the  p-region. 

/s 

For  the  density  operator  r^’(f),  separating  the  effects  of  the  external  light  field  from 
the  internal  coupling,  as  before  we  have 


+ (13) 

where  rQ(f)  represents  the  p-region  density  operator  in  the  absence  of  the  light  field,  and 
rj  {t)  represents  the  p-region  density  operator  changes  due  to  the  light  field  interaction 
with  the  extended  system  and  the  consequent  changes  in  electronic  density  distribution. 
The  effective  Hamiltonian  for  the  s-region  now  can  be  written  as, 

F.{t)  = H,  + H[(t)  + &i(t)  + G'jt) 

G2(()  = (rp[fg(f)ff,,,(f)]/(r„[f'>(()]  (14) 

where  comes  from  the  coupling  between  the  s-region  and  the  p-region  in  the 

absence  of  the  light  field  and  G[^{t)  comes  from  the  effects  on  the  s-region  of  the 
light  field  interactions  with  the  extended  system.  The  coupled  SCF  equations  could 
be  constructed  if  complete  information  were  known  about  the  s-region.  However  the 
s-region  is  described  only  through  its  statistical  averages  in  this  treatment.  To  introduce 
these  averages,  we  rearrange  the  equation  for  the  p-region  to  make  explicit  the  fluctuation 
and  dissipation  forces,  and  interaction  terms  between  the  light  field  and  the  extended 
system  which  depend  only  on  the  p-region  variable.  We  derive  this  effective  Hamiltonian 
of  the  p-region  in  the  next  section. 
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Self  Consistent  Field  (SCF)  Treatment  of  the  Coupling  Regions 


The  Liouville-von  Neumann  equation  for  the  p-region  is 

dt‘Ut) 


dt 


= [H,  + H'M  + G»(()  + G'M,  ri’(()l 


P 


P 


= Fwi  + Iffiif).  f'’(«)l  + [gS((),  f'’(()]  + [GL(t). 


P 


P 


If  we  define, 


Hp  - [Hp , •] 


s°{t)  = ig;(i),h 


P 


we  can  formally  write  the  Liouville-von  Neumann  equation  for  the  p-region  as, 


(15) 


(16) 


i — r — = [iip  + Hp{t)  -t-  G„{t)  + Q {t)]Y^{t) 


dt 


(17) 


Similarly,  the  Liouville-von  Neumann  equation  for  the  s-region  is,  with  corre- 
sponding definitions. 


— [Hfi  + '^i(^)  + ^s(0  + ^s(^)]f'‘’(0 


(18) 


We  have  to  solve  Liouville-von  Neumann  equations  (17)  and  (18)  simultaneously. 
Here  we  assume  we  already  know  the  density  operator  for  the  s-region,  so  we  can 
proceed  to  find  the  effective  Hamiltonian  of  the  p-region.  In  the  interaction  picture  we 
can  rewrite  equation  (17)  as 


= w;(()7f'’(J)/ 


(19) 


with 


= exp  (^-iipt^[fi'p{t)  -f-  Qp{t.)  -t-  Ql{t)]  exp  [-ilipt 


(20) 


17 


and 


(21) 


The  equation  (19)  can  be  formally  integrated  over  time,  and  after  transforming  back 
to  the  Schroedinger  picture,  the  p-region  density  operator  is  given  by 


r^{t)  =exp  ( )f^(o) 


t 

i I exp  - /'))[^;,(^')  + g^p{f')  + gl{t')]tp{t')df' 

b 


(22) 


Now  if  we  put  the  solution  Eq.  (22)  back  into  the  last  term  on  the  right  hand  side 
of  equation  (17),  the  result  is  an  integrodifferential  equation  for  the  p-region  density 
operator. 


d 


= n,P{t)  + nl{t)rp{t)  + (of^(f) 


-\-Qp{t)  exp(-?;?fpt)r(o) 


(23) 


I 


^0/^/ 


0 


After  a few  simplifications,  expression  (23)  can  be  written  as  a generalized  Langevin 
expression  (GLE)  with  three  terms  due  to  the  interaction  of  light  with  both  the  s-  and 
p-  regions, 


0 

i-^v\t)  = H,P{t)  + n[{t)t^{t)  + gl{t)rP{t)  + np{t)rp{o) 

f f 

-i  I Mp{t,t')P{t')d£ -i  f M^j,{tj')FP{t')d.t' 

0 b 


(24) 


The  first  term  on  the  right  hand  side  refers  to  the  Liouville— von  Neumann  equation 
of  the  p-region  in  the  absence  of  the  s-region  and  the  light  field.  The  second  term  gives 
the  interaction  between  the  light  field  and  the  p-region  in  the  absence  of  the  s-region. 
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The  third  term  gives  the  effective  effects  on  the  p-region  due  to  the  interaction  of  the 
light  field  with  the  s-region.  The  fourth  term  can  be  written  as 


(25) 

and  is  the  “energy  fluctuation”  term  as  in  the  GLE  equation  without  the  light  field. 
The  “memory  kernel”  of  the  fifth  term  given  by 

= Ql{t)  ex\')[-iiip{t  - t')]Qp{t')  (26) 

is  the  “energy  dissipation”  term  in  the  GLE  equation  without  the  light  field.  It  is 
obvious  that 

Mp{tj')  =np{t)^{t')  (27) 

The  “memory  kernel”  of  the  sixth  term  for  coupling  between  the  internal  coupling 
term  and  the  light  field  coupling  term  is  given  by 


The  density  operator  T^(f)  depends  on  the  temperature  and  the  initial  values  of 
properties  in  the  s-region.  Treating  the  s-region  as  a stochastic  environment,  we  can 
take  the  average  of  the  previous  equation  over  initial  conditions  using  r‘’(0)  and  assume 


that  [43]: 

(i)  the  fluctuation  force  averages  to  zero  on  the  primary  time  scale,  i.e. 


{Upit))  = 0 


(29) 


(ii)  the  memory  kernel  describes  both  instantaneous  and  delayed  dissipation  leading 
to  a temporary  equilibrium,  so  that 


1 


= -Wp{t)S[t  - t')  +J^p(t)7][t  - t') 


(30) 
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where  the  first  term  to  the  right  contains  a Dirac  delta  function  and  describes  instanta- 

/S 

neous  dissipation  while  in  the  second  term  the  operator  Mp  describes  a delayed  dissi- 
pation with  a time  delay  implicit  in  the  analytical  function  r/. 

This  expression  for  the  memory  kernel  is  general  enough  to  describe  most  cases  of 
interest.  In  the  present  study  we  are  also  assuming,  and  supported  by  experimental  evi- 
dence, that  the  response  of  the  s-region  electronic  system  is  very  fast,  and  consequently 
we  shall  keep  only  the  first,  instantaneous  dissipation  term  from  now  on.  Dissipative  po- 
tentials  Wp  can  be  calculated  within  a model  or  parametrized  to  reproduce  experimental 
desorption  rates. 

Adding  the  third  assumption  that  the  factors  in  the  memory  kernel  with  and  without 
the  light  field  are  large  at  different  times,  we  also  write 


(31) 


Using  the  three  basic  assumptions  above,  we  can  rewrite  Eq.  (24)  as 


✓N 


(32) 


with 


VVp(i)  = .] 


(33) 


Returning  to  the  Liouville-von  Neumann  equation  for  the  p-region,  we  have 


(34) 


with  the  effective  Hamiltonian  of  the  p-region  given  by 


/S  /s 


Fp{t)  = Hp  + Hl{t)  + G’Jt)  - iWp{t)/2 


P 


P 


(35) 
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The  primary  Hamiltonian  Hp,  which  is  the  sum  of  the  kinetic  energy  and  the  potential 
energy  operator  and  is  dependent  on  the  p-region  degree  of  freedom  with, 

Hp  = kp  + Vp  (36) 

and  the  direct  light  interaction  term  Hp{t),  can  be  obtained  directly  from  electronic 
structure  calculations  of  the  p-region. 

The  indirect  light  interaction  term  and  the  dissipation  term  must  be  derived  from 
the  explicit  expressions  of  the  secondary  region.  A detailed  treatment  of  the  s-region 
with  driving  forces  from  p-region  and  the  external  light  field  will  be  given  in  Ch.  4. 

Direct  and  Indirect  Photodesorption 

Depending  on  whether  the  laser  pulse  photons  are  absorbed  by  the  adsorbate 
molecules,  or  absorbed  by  the  substrate,  the  photodesorption  processes  can  be  catalogued 
into  two  types:  direct  and  indirect  electronic  excitation  photodesorption. 

In  direct  photodesorption,  there  are  several  possible  mechanisms  by  which  absorbed 
molecules  can  be  excited  [50]:  by  resonant  valence  excitation  without  changing  the 
molecule’s  charge  states  [51];  through  photoionization,  resulting  in  a positively  charged 
excited  state,  as  proposed  for  CO/Pt(lll)  [52,  53];  theoretically,  if  the  hole  created 
by  direct  valence  excitation  of  an  absorbed  molecule  is  filled  by  a substrate  electron, 
the  excited  state  can  also  correspond  to  a negatively  charged  state;  another  possible 
direct  excitation  mechanism  is  the  photoexcitation  of  the  adsorbate-substrate  bond,  e.g., 
transition  of  an  electron  from  a bonding  to  an  antibonding  orbital.  Such  an  excitation 
process  may  or  may  not  have  charge  transfer  character. 

Here  we  briefly  describe  the  direct  mechanism  we  used  in  previous  calculations  by 
Yi  and  Micha[44]  for  CO  photodesorption  from  Ni(OOl). 
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The  effective  Hamiltonian  for  the  direct  interaction  of  the  light  field  with  the  p- 
region  is  given  by 


Fp{t)  = Hj,  + Hlit)  - iW,{t)/2 


(37) 


By  using  the  dipole  approximation  for  the  interaction  term,  we  obtain 


?, 


F,  = ffp  - S{t)  ■ D,  - -W„ 


(38) 


where  £{t)  is  the  electric  field  and  Dp  is  the  transition  dipole  in  the  p-region. 

Using  time-dependent  perturbation  theory  to  treat  both  the  dissipation  and  light 


field  with 


Fp  = f; + F,:(() 


F„‘(t)  = -S(t)  ■ Dp  - -W. 


I 


(39) 


p 


p 


separated  into  zeroth-  and  first-order  terms,  and  working  with  matrices  in  the  diabatic 


electronic  basis,  for  the  special  case  of  a initial  pure  state  |0o(^))>  density  operator 
is  given  by 


f'-(t)  = mmrm 


(40) 


So  the  resulting  state  is  given  to  first  order  by 


with 


'ipP{t))  = [ijjP[t))+i  I ds  f(.s)exp[-?;F^^s]Dp|-(/>g(.s)) 

b 


+i  I ds  exp[-iFjjs]^Wp\ij^{s)) 

b 


(41) 


•0g(s))  = e,'rp[-?;F°.s](0g(O)) 


(42) 


keeping  in  mind  that  ipp  is  not  an  eigenfunction  of  Fjj,  and  is  obtained  by  constructing 
the  time  integrals  while  propagating  the  dipole-driven  and  dissipation-driven  transition 
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amplitudes  from  the  first  and  second  integrals  respectively.  It  should  be  noted  that  we 
can  solve  the  direct  photodesorption  problem  numerically  without  using  perturbation 
theory. 

In  indirect  photodesorption,  electronic  excitations  created  in  the  substrate  by  photon 
absorption  travel  to  the  surface  and  induce  reaction  by  localizing  on  the  adsorbate  or 
the  adsorbate-substrate  complex.  In  most  cases,  photon  absorption  creates  energetic 
electron-hole  pairs  through  bandgap  excitation  as  the  initial  electronic  excitations. 
Photodesorptions  from  metal  surfaces  have  been  attributed  to  photoexcited  or  “hot 
electrons”,  such  as  NO  photodesorption  from  Pd(lll)  [28,  30],  CO  photodesorption 
from  Cu(lll)  and  Cu(OOl)  [29,  39,  54],  and  O2  photodesorption  from  Pt(lll)  [34]. 
From  now  on,  we  will  focus  on  indirect  photodesorption  because  it  predominates  in 
most  photodesorption  processes  [25].  By  studying  the  interactions  between  the  light  field 
and  the  substrate  metal  surface  and  its  coupling  with  the  p-region,  we  will  introduce 
an  effective  transition  dipole  for  the  p-region  which  is  explicitly  dependent  on  the 
light  pulse  fluence.  Based  on  the  effective  transition  dipole,  we  will  develop  a new 
indirect  photodesorption  mechanism  and  apply  it  to  the  CO/Cu(001)  photodesorption 
process.  Formally,  for  indirect  photodesorption  in  the  dipole  approximation,  the  effective 
Hamiltonian  for  the  p-region  is  given  by 

F,{t)  = Hj,  - S{t)  ■ ip{t)  - iWp{t)/2  (43) 

/s 

where  is  the  effective  transition  dipole  for  the  p-region. 


CHAPTER  3 

DYNAMICS  OF  THE  PRIMARY  REGION 
Time-dependent  Schroedinger  Equation  for  the  Primary  Region 

In  this  chapter,  we  study  the  p-region  dynamics  of  the  extended  system.  The 
effective  Hamiltonian  of  the  p-region  , derived  from  the  SCF  treatment  of  last  chapter 
is  given  by, 

Fp{t)  = Hp  + Hl(t)  + G\,{t)  - iWp{t)/2  (44) 

with  the  primary  region  density  operator  written  given  by 

(45) 

i 

where  is  the  weight  function  of  each  pure  state.  We  can  write  the  time  dependent 
Schroedinger  equation  (TDSE)  of  the  primary  region  for  each  pure  state  as  (drop  the 
index  i) 

= Fp{t)\ipP{t))  (46) 

where  h = I here  and  henceforth.  The  wavefunction  4>P(t)  is  a function  of  all  the 
p-region  coordinates  including  nuclei  and  electrons.  It  fully  describes  the  dynamics  of 
the  adsorbate  in  the  presence  of  the  substrate  and  the  light  field.  Instead  of  trying  to 
solve  the  time  dependent  Schroedinger  equation  for  all  the  particles  simultaneously,  we 
adopt  the  Born-Oppenheimer  approximation,  which  takes  note  of  the  great  difference  in 
masses  of  electrons  and  nuclei.  Because  of  this  difference,  the  electrons  can  respond 
almost  instantaneously  to  displacement  of  the  nuclei.  Therefore,  it  is  possible  to  regard 
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the  nuclei  as  fixed  in  position  and  to  solve  the  TDSE  for  the  electrons  in  the  static  electric 

potential  arising  from  the  nuclei  in  that  particular  arrangement.  Different  arrangements 

of  nuclei  may  then  be  adopted  and  the  calculation  repeated.  The  set  of  solutions  so 

obtained  allows  us  to  construct  the  potential  energy  surface  (PES).  We  can  expand  the 

p-region  wave  function  in  a basis  of  electronic  states  {|0j(Q))}  which  parametrically 

depend  on  nuclear  variables  in  the  p-region.  The  coefficients  of  the  expansion  are  then 

nuclear  wave  functions  which  only  depend  on  nuclear  variables  and  time, 

N 

|■0'’(Q,'))  = ^«(Q.^)IV>5(Q)>  (47) 

J=1 

where  Q = (Qi,  Q2,Qi,  ■ ■ ■)  is  a set  of  generalized  coordinates  for  the  motions  of  the 
nuclei.  The  p-region  effective  Hamiltonian  is  given  by 

I,j  (48) 

^ /s  '^7  '^7  7t 

where  K[  j,  V[j,  Hj  j,  Gj  j and  Wjj  are  operators  only  depending  on  the  p-region 
nuclear  variables 

The  expansion  over  this  basis  is  then  introduced  into  the  TDSE,  so  this  can  be 
written  in  the  matrix  form 


r\ 


(49) 


with 

F"(Q4)  = K+V(Q)+H'(Q,  f)+G'(Q,  t)-'-W{Q,  t)  (50) 

where  K is  the  kinetic  energy  dependent  on  nuclear  variables,  a N by  N matrix,  V,  H^, 
G and  W are  the  N by  N matrices  of  the  complex  potential  energy  terms  dependent  on 
nuclear  variables,  and  1/7"  (Q,/)  is  a nuclear  wavefunction,  a N by  1 column  matrix. 
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The  elements  of  the  kinetic  energy  matrix  K are 


1 


2M 


T 2 


OQ 


(9Q 


(51) 


where  the  second  term  in  the  square  bracket  represents  the  momentum  coupling  between 
the  electronic  basis  states. 

The  elements  of  the  real  potential  energy  matrix  V are  given  by 


[v(Q)l;,j  = W(Q)|l>{Q)|V’5(Q)> 


(52) 


and  they  represent  all  Coulomb-type  interactions  between  the  electronic  basis  states 
The  elements  of  the  real  potential  energy  matrix  are  given  by 

(H'{Q,<)I/,J  = {</’f(Q)|ffJ(Q.*)IV'5(Q)>  (53) 

and  they  represent  the  direct  interactions  between  the  electronic  basis  state  and  light  field. 
The  elements  of  the  real  potential  energy  matrix  are  given  by 

[G'(Q, <)]/,/  = {V’f(Q)|G!,(Q.')l'/'j{Q)>  (54) 

and  they  represent  the  effective  interactions  between  the  electronic  basis  state  and  light 
field  which  is  from  the  light  field  interaction  with  the  s-region. 

Finally,  the  imaginary  potential  matrix  elements  are  given  by 

(W(Q,t)l/,j  = (■0/(Q)Fp(Q.<)IV'5(Q))  (55) 

We  can  develop  our  treatment  of  the  molecular  dynamics  in  any  convenient  elec- 
tronic representation.  Two  of  them  are  frequently  used:  the  adiabatic  (a)  representation 
where  V is  diagonal  [55-57],  and  the  diabatic  (d)  representation  where  K is  diagonal. 
Electronic  structure  calculations  provide  results  in  the  adiabatic  representation,  but  in- 
tegration of  the  dynamics  equations  is  simpler  in  the  diabatic  representation.  Therefore 
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we  work  with  matrices  in  a diabatic  basis,  constructed  so  that  transforming  back  to  the 
adiabatic  basis  recovers  the  calculated  electronic  potential  surfaces. 

We  developed  the  physical  model  in  the  diabatic  representation  which  means  there 
are  no  coupling  terms  in  kinetic  energy  matrix.  The  elements  of  all  matrices  are  all 
functions  of  the  primary  coordinates  and  time,  and  all  of  them  are  real  and  symmetric 
matrices.  This  means  that  they  can  be  diagonalized  by  orthogonal  transformations.  For 
the  real  part  of  the  potential,  the  diagonalization  of  is  in  fact  the  transformation 
into  the  adiabatic  representation,  as  defined  for  real  potentials[56,  57].  It  is  given  by 


V(^>(Q)=LvV(‘i>(Q).LV 


(56) 


where  the  orthogonal  transformation  matrix  Ly  is  an  element  of  the  SO(N)  algebra  [58]. 

In  this  study,  we  assume  the  electrons  can  be  only  on  two  states;  ground  and  excited 
electronic  states.  The  matrices  in  the  electronically  diabatic  representation  are 


k('^)  = L_^i 


Gf{Q,t)  = 


2MdQ‘ 

V,{Q)  Vc{Q) 

LK-(Q)  K(Q)J 

HiiQjy 

mQ,t)  G[{Qj.y 
GiiQJ.)  G[{Q,t)_ 

'WgiQj.)  WciQJ,) 
Wc(Q,t)  We{Q,t) 


(57) 


(58) 


(59) 


(60) 


(61) 


In  these  expressions,  I is  the  identity  matrix,  g represents  the  electron  diabatic  ground 


state,  e represents  the  electron  diabatic  excited  state  and  c represents  the  coupling 


between  these  two  states. 
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We  now  give  the  transformation  matrix  between  adiabatic  and  diabatic  representa- 
tions for  the  two  electronic  level  system,  the  transformation  matrix  Ly  can  be  written  as 


cos7;(Q)  siii7/(Q) 
— sinr;(Q)  cos7/(Q) 


with  a rotation  angle  that  is  a function  of  primary  coordinates, 


viQ) 


1 

- arctan 
2 


2K.(Q) 

Ve{Q)-V,{Q) 


We  treat  the  interaction  between  the  light  field  and  the  system  within  the  electric 
dipole  approximation,  for  direct  photodesorption,  with  a semiclassical  treatment  for  the 
light  field,  so  that  we  have 


= D{Q)-£{t) 


(64) 


The  matrix  for  the  electric  dipole  operator,  needed  to  describe  transitions  induced  by 
the  direct  light  field  effect,  is  given  by 


D('^>(Q) 


'D^(Q)  DciQ) 
Dc{Q)  De{Q) 


(65) 


In  contrast  to  the  direct  photodesorption,  for  indirect  photodesorption,  we  have 


= A {Qj,)-£{t) 


(66) 


The  matrix  of  the  effective  transition  dipole  which  depend  only  on  p-region  nuclear 
variables  and  light  field  fluence  is  given  by 


A('^((Q,f) 


A^{Q,t)  Ac(Q,0 

AciQJ,)  Ae{Q,t) 


The  effective  transition  dipole  comes  from  the  interactions  between  the  light  field  and 
s-region  electron-hole  excitations.  Note  that  it  depend  on  p-region  nuclear  variables 


28 


and  the  time  variable.  In  the  next  chapter,  we  will  give  a detailed  discussion  about  this 
effective  dipole  operator. 

We  now  give  a brief  discussion  on  how  to  get  those  matrices  in  the  electronic 
diabatic  representation.  A more  detailed  discussion  will  be  given  in  Ch.  4 for  the 
dissipation  term  and  the  effective  transition  dipole  term  of  the  p-region,  and  in  Ch.  6 
for  the  case  of  CO  photodesorption  from  Cu(OOl). 

For  matrix  elements  of  V^"^(Q),  we  can  perform  semiempirical  molecular  orbital 
calculations  using  parameterization  of  electron  integrals  with  incomplete  neglect  of  dif- 
ferential overlap  appropriate  for  systems  containing  transition  metal  ions  (the  ZINDO 
parametrization)  [46-48].  The  ground  electronic  state  was  obtained  from  ZINDO/ 1 cal- 
culations, which  are  parametrized  to  provide  the  ground  electronic  state.  Electronically 
excited  states  were  obtained  from  ZINDO/S  calculations.  These  construct  many-electron 
states  as  linear  combinations  of  singly  excited  Slater  determinants  obtained  from  a ref- 
erence determinant. 

For  the  matrix  elements  of  the  primary  dipole  D*''^(Q),  we  know  that  Dg=Dg=0 
because  of  selection  rules  on  parity.  When  we  get  the  electronic  excited  states  in 
ZINDO/S  calculations,  they  also  gives  us  the  transition  dipoles  between  the  excited 
states  and  ground  state  simultaneously. 

For  the  dissipation  matrix  elements  of  W^'^*(Q,/),  we  shall  formally  give  in  Ch.  4 
the  explicit  form  of  the  dissipation  dependence  on  properties  of  the  p-  and  s-  regions. 

Finally,  we  shall  obtain  the  effective  transition  dipole  from  both  linear 

response  and  nonlinear  response  expressions,  treating  the  metal  substrate  as  a continuous 
medium  within  first  order  perturbation  theory  and  including  all  order  response. 


29 


Physical  Properties  of  Interest 

As  the  wavefunctions  are  calculated  by  propagation  over  time,  it  is  possible  to  use 
them  to  calculate  time  dependent  properties.  Here  we  give  the  expressions  which  we 
shall  use  in  Ch.  6 and  Ch.7  for  the  CO/Cu(001)  system.  Our  model  deals  with  a CO 
molecule  adsorbed  with  its  axis  perpendicular  to  the  surface,  on  top  of  a Cu  atom, 
with  the  coordinate  Z representing  the  distance  from  the  center  of  mass  of  CO  to  the 
Cu  surface.  The  dynamics  of  desorption  was  assumed  to  involve  the  softest  vibrational 
mode  of  the  absorbed  CO,  which  is  its  frustrated  translation  mode,  or  “T  mode”,  parallel 
to  the  surface  with  coordinate  x representing  displacements  from  the  vertical  axis. 

We  consider  the  experimental  situation  where  the  system  is  initially  in  its  ground 
electronic  state  g,  in  the  ground  vibrational  state  in  the  Z motion  and  a certain  vibration 
state  Vx  in  the  x motion,  from  which  it  is  excited  by  a light  pulse. 

In  prevous  works,  for  the  CO/Ni(001)  system  [44],  assuming  a direct  photodes- 
orption mechanism,  we  calculated  the  time  evolution  of  vibrational  population  distri- 
butions on  both  ground  and  excited  electronic  potential  energy  surfaces  [44].  For  the 
CO/Cu(001)  system,  assuming  a indirect  photodesorption  mechanism  and  using  pertur- 
bation theory,  we  calculated  the  time  evolution  of  survival  probabilities  with  several 
dissipation  strengths  [45]. 

In  this  study,  we  shall  concentrate  on  the  time  dependent  population  distributions 
on  different  potential  surfaces  and  the  time  dependent  desorption  yield.  We  shall  give 
detailed  results  for  their  dependence  on  the  light  pulse  strength,  shape,  the  substrate 
dissipation  and  initial  temperature,  and  the  adsorbate  isotope  masses. 

The  desorption  yield  can  be  obtained  from 

= lim  Yr^{t) 

f-^OC 


(68) 
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where 


CX3  ^ max 


i.e.  by  integrating  the  state  density  on  both  electronic  ground  and  excited  potential 
surfaces  from  a desorption  distance  to  infinity  and  including  all  the  x component, 
and  taking  the  infinite  time  limit.  If  we  consider  the  initial  temperature  of  the  substrate, 
the  desorption  will  depend  on  the  initial  thermal  vibrational  state  and  the  desorption 
yield  Y will  actually  be  given  by 


y = Yl  (^0) 

where  is  a state-resolved  desorption  yield  obtain  from  Eqs.  (68)  and  (69)  and 
is  a weight  function.  If  the  substrate  surface  is  at  temperature  T,  before  a light  pulse 
arrives,  is  the  Boltzmann  function 


exp{-eyJkBTs) 

'Z^expi-eyJk.DT,) 


(71) 


Here,  €y^  is  the  vibrational  energy  of  level  Vj-,  and  ku  is  the  Boltzmann  constant.  From 
this  expression  we  can  study  how  the  initial  substrate  temperature  affects  the  desorption 
yield. 

The  time  dependent  populations  on  the  electronic  ground  and  excited  potential 
surfaces  are  given  by 


^T,vA^)=  (I  = e,g) 


(72) 


Similarly  to  the  desorption  yield  calculations,  if  we  consider  the  initial  temperature  of 
the  substrate,  the  vibrational  state  population  will  depend  on  the  temperature  and  the 
population  distribution  Pi{t)  will  actually  be  given  by 

P,{t)  = 

Vx 


(73) 
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From  the  calculations  of  this  property,  we  shall  see  that  there  are  multiple  electronic 
excitations  and  deexcitations  during  the  desorption  process. 


CHAPTER  4 

GENERAL  THEORY  OF  THE  SECONDARY  REGION 

Introduction 


This  chapter  presents  the  general  theory  of  the  s-region  and  its  couplings  with 
both  the  p-region  and  the  time  dependent  light  field.  The  s-region  can  respond 
to  these  driving  forces  with  lattice  vibrations  and  electronic  excitations.  We  note 
that  the  lattice  vibrations  and  electronic  excitations  have  very  different  time  scales; 
the  lattice  vibration  excitations  have  time  scales  on  the  order  of  nanoseconds  while 
the  electronic  excitations  have  time  scales  in  the  femtosecond  range.  Because  we 
are  interested  in  the  femtosecond  dynamics  in  the  s-region,  we  ignore  the  substrate 
vibration  excitation  effects  in  the  following  discussion.  The  photodesorption  mechanism 
including  the  resonance  absorption  of  light  by  the  substrate  metal  surface  leading  to 
electronic  excitation,  the  transfer  of  electronic  energy  into  the  adsorbate  and  electronic 
rearrangement  in  the  primary  region  coupled  to  dissipation  into  the  secondary  region 
is  illustrated  in  Fig.  (3) 

The  effective  Hamiltonian  for  the  s-region  from  Ch.  2,  when  we  ignore  the 
interaction  between  the  light  field  and  the  p-  and  s-  region  coupling  term,  is  given  by 


H.  + + G<>M 


We  recall  that  all  the  terms  in  Eq.  (74)  include  only  charges  in  the  s-region.  In  section 
1 we  will  discuss  the  coupling  force  effects  on  the  s-region  and  give  the  dissipation 
term  for  the  p-region.  In  section  2 we  will  discuss  the  interaction  of  the  light  field  with 
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the  s-region  using  perturbation  theory  and  give  its  linear  and  nonlinear  response  and  the 
effective  transition  dipole  for  the  p-region. 


Dissipation  Effect  of  the  Secondary  Region 


The  calculation  of  the  dissipation  potentials  of  the  p-region  requires  the  calculation 
of  electron-hole  excitations  in  the  s-region.  This  can  be  done  with  a large  cluster  of 
metal  atoms,  and  has  been  carried  out  for  the  ground  electronic  state  [20-22].  In  this 
study  we  need  the  dissipation  potentials  in  the  diabatic  representation  for  the  ground  and 
excited  states  and  their  coupling.  An  expression  has  been  derived  for  them  by  Beksic 
and  Micha  [43].  In  this  section,  we  derive  expressions  for  the  dissipation  potentials 
using  the  density  matrix  approach. 

The  effective  Hamiltonian  for  the  electron-hole  excitations  created  in  the  s-region, 
when  this  is  driven  by  forces  from  the  p-region,  can  be  written  using  creation  and 
annihilation  operators  as 


F,{t)  = H,  + G%.) 


A A A' 

where  the  first  summation  describes  the  electrons  in  states  A = {h,k)  for  the  energy 
band  b,  and  given  electron  wavevector  k.  The  double  summation  describes  electronic 
transitions  induced  by  external  forces,  with  the  factor  representing  an  energy 

change  due  to  electron  scattering.  It  should  be  noted  that  the  driving  force  due  to 
the  coupling  between  the  p-region  and  s-region  can  have  a wide  range  of  energies. 
The  transitions  occurring  in  the  s-region  can  be  either  intraband  or  interband  transitions. 
Interband  transitions  can  be  classified  as  phonon-assisted  or  vertical  interband  transitions 
corresponding  to  Ak  ^ 0 and  AA:  = 0 respectively. 
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The  p-  and  s-region  coupling  Hamiltonian  for  the  electron  can  be  written  as 


✓N 


where  we  have  expanded  it  in  the  electron  operators  and  the  coefficients  Hly{t)  are 
yet  operators  in  the  p-region.  This  coupling  therefore  contains  long  range  Coulomb 
interactions  between  charges  of  the  adsorbate  and  the  substrate,  and  diabatic  couplings 
from  the  elimination  of  momentum  couplings. 

From  Ch.  2,  we  know  that 


(77) 


If  we  substitute  Eq.  (76)  into  the  above  equation,  we  obtain 


I /N 


A A' 


Comparing  with  Eq.  (75),  we  see  that 


i7A,A'(^) 


fr,,[fp(/)] 


(78) 


(V9) 


If  only  one-electron  processes  are  of  interest,  the  N-electron  density  can  be  reduced  to 


the  one-electron  density  matrix,  which  is  defined  as 


ff{t)  = fr2...A-[P(f)]  = '^'^px,y{t)axa\, 

A A' 


(80) 


If  we  put  Eqs.  (75)  and  (80)  into  the  L-vN  equation  for  the  s-region,  we  obtain 


The  one-electron  density  matrix  elements  pa,A'(0  ^^e  evaluated  in  the  basis  of  electronic 
states  I A)  of  the  isolated  solid  without  the  adsorbate.  After  evaluating  the  commutator 
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on  the  right-hand  side  of  Eq.  (81),  With  some  simplifications,  we  obtain  the  following 
first  order  differential  equation  in  time  that  couples  the  density  matrix  elements. 


dt 


+ — PA,A"(f)l7A",A'{^)] 


A 


// 


(82) 


These  equations  of  motion  can  be  solved  using  the  metal  dielectric  function  if  the 
coupling  is  weak  [59].  If  the  coupling  is  strong,  it  can  be  solved  by  temporal  linearization 
[60].  In  the  zeroth  order  approximation,  which  means  no  coupling  with  the  p-region, 
the  solution  is  oscillatory  and  is  given  by 


pIa'(^)  = exp[-7(eA  - £a')^]Pa,A'((^) 


0 


(83) 


Since  the  s-region  is  described  by  one  electron  density  operator,  we  can  now  replace 
the  s-region  reduced  density  operator  into  Gp{t)  from  Ch.  2,  to  obtain 


tr,,[t^Q{t)Hps{t)] 


(84) 


A A' 


where 


7a,A'(^)  = 


PA,A'(^) 

Epa,a(0 

A 


(85) 


Replacing  this  Gp{t)  into  the  energy  fluctuation  and  memory  terms  from  Ch.  2 
and  performing  the  statistical  average  ( ) over  the  distribution  of  initial  conditions  in 
the  s-region,  we  obtain 


{Tipit))  = {g^{t))  exp(-mpt) 


p 


(7A,A'(f))^A,A'(^)  exp(-74fpf) 


(86) 
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and 


A, A'  A", A'" 


(87) 


in  the  notation  of  Liouville  operators. 

These  expressions  are  similar  to  the  ones  in  the  literature  [61]  and  can  be  analyzed 
as  Brownian  motion.  The  relative  density  amplitudes  7a,a'(^^)  are  functions  of  the  set  of 
their  initial  values  {7a,a'(0)}  with  a probability  distribution  over  which  the  fluctuation 
and  dissipation  terms  must  be  averaged.  From  Eq.  (82)  we  know  that  the  7a,A'(0 
oscillate  rapidly  and  irregularly  on  the  time  scale  of  the  adsorbate  motions.  Therefore 
we  assume  that  for  all  times  ^ > 0, 


(7a,A'(^))=0  (88) 

and 

(7a,A'(^)7A",A'"(J^  + ■’■))  = 27T()'(r)CA,A',A",A"'(^)  (89) 

where  r — t'  — t and  the  r-dependence  on  the  right  indicates  that  the  instantaneous 
correlation  involves  a temporary  equilibrium  distribution. 

We  find,  integrating  over  r both  sides,  that 

OC 

C'A,A',A»,A'''(f)  = ^ j dT{'yx,X'{thx",\>"{t  + T))  (90) 

— C5C 

The  fluctuation  and  memory  terms  can  now  be  written  as 

{n{t))  = o (91) 

and 


{Mp{t,t'))  = 6{T)Wj,{t)/2 


(92) 
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with 


W,(t}  = inY,  E Ca, A', 

A, A'  A", A'" 


(93) 


This  is  similar  to  the  expression  given  by  Beksic  and  Micha  [43]  but  expressed  in  terms 
of  Liouville  operators.  It  shows  how  dissipation  related  to  change  density  fluctuations 
in  the  s-region.  A corresponding  expression  can  be  derived  working  with  operators 
in  Hilbert  space  instead  of  superoperators  in  Liouville  space,  and  introducing  the 
time-dependent  equations  of  motion  for  the  amplitudes  0^’.  This  leads  to  a dissipation 


operator 


Wp{t)  = 47t  EE  CA,A.,A”A”(<)Ht,j,(i)«j;.,A".(i) 

A, A'  A", A'" 


(94) 


If  we  expand  Eq.  (94)  in  the  p-region  diabatic  electronic  states  introduced  in  Ch. 
3,  we  obtain  the  expression 

A, A', A", A"'  K ’ ’ (95) 


where  we  find  the  time-correlation  function  of  the  electronic  density  fluctuations  and  the 
product  of  Hamiltonian  matrix  elements  {Q,  t)  of  Coulomb  interactions  between 
p-  and  s-regions,  which  involve  transitions  [KX)  {IX')  . These  therefore  contain 
the  overlaps  of  p-  and  s-  region  electronic  states  and  are  largest  when  I =K.  These 
functions  of  Q decrease  at  large  distances  and  can  be  written  as  sums  of  short  and  long 
range  terms.  Here  we  shall  not  elaborate  on  them,  and  we  will  introduce  parametrized 
dissipation  functions  in  the  calculations. 
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Interaction  of  a Light  Pulse  with  the  Secondary 
Region;  Linear  and  Non-linear  Response  Theory 

In  this  section,  treating  the  metal  surface  as  a continuous  medium,  we  picture  the 
electric  field  penetrating  the  metal  substrate.  We  know  that  it  undergoes  dispersion  and 
absorption  so  that  its  amplitude  decreases  exponentially  with  increasing  depth  2.,  inside 
the  metal.  A local  field  parallel  to  the  surface  with  absorption  and  dispersion  can  be 
expressed  as 


OC 

./  (96) 

— OC 

£loc-A^s)  - S.^exp[-'^{u})z,^]exp[iK{uj)z.,,] 

where  k and  7 are  the  dispersive  wavevector  and  absorption  coefficients  respectively. 
This  field  drives  electronic  excitations  inside  the  metal  by  coupling  to  the  dipole  density 
arising  from  electronic  fluctuations  in  the  s-region.  The  interaction  between  light  pulse 
and  substrate  metal  surface  is  illustrated  in  Fig.  (4) 

Positioning  the  semi-infinite  metal  with  its  surface  perpendicular  to  the  z-axis  at 
Zs  = 0 and  extending  to  negative  values,  and  employing  a secondary  dipole  per 
unit  length  £>,(2,),  we  write  the  effective  Hamiltonian  in  the  s-region  in  the  dipole 
approximation  as 


F.(>) 


0 

I dz  Df,{z,o)£i„c{zs,t) 

— OC 


From  Ch.  2,  we  know  that  the  density  operator  of  the  s-region  can  be  separated  into 
two  parts 


r*(()  = rs(«)  + r?(i) 


with  the  initial  conditions  at  time  t = 0,  f ‘'(0)  = Tq. 
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Figure  4;  Interaction  between  light  pulse  and  the  substrate  metal  surface 
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In  the  interaction  picture,  we  can  write  the  L-vN  equation  as, 


orHt)T 


(99) 


with 


where  the  zeroth-order  evolution  operator  is 


(100) 


U^{t)  = exp{—iHf,t) 


(101) 


We  can  formally  solve  Eq.  (99)  in  the  interaction  picture  and  transfer  it  back  to  the 
Schroedinger  picture;  we  obtain  to  all  orders  in  the  perturbing  field 


r-'(f)  = r5(()  - i I dt' 

b 

The  perturbation  term  for  a weak  field,  can  be  written  to  the  first  order  as 


= -i  I dt'  (fht)[upt'yHi(t')opt'),n]u«(t) 


/S  /S 


0/^xt 


= 'I 


— oc 


0 

0 f 

b 


dz. 


t 


(102) 


(103) 


Replacing  the  perturbed  s-region  density  operator  in  the  p-region  effective  Hamiltonian, 


Fp{t)  = H,  + G'^{t) 


(104) 


with 


G^it)  = «r,[f;(()//p,l/«r,[f-'(()I 


(105) 
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we  obtain 


0 


G'(t)  = 


(r,{[  / dz,  fdf  UO{t)lU»(tn^D,(z.)U!!(t'}£,„{z.,t),  f5](7?(t)'|ff,,,} 

0 


— oc 


0 


- / r/2s  / dt'Ap^\tj':z^)£j„c{z<,j') 


— OC 


0 


(106) 


with 


2s)  = i 


0/4\t 


(107) 


/N 

where  Ap  is  the  effective  dipole  density  created  in  the  p-region  by  its  indirect  response 
to  the  field  in  the  substrate,  and  is  given  explicitly  by 


D,{z,j')j  = U^(t'yD,{z,)U^{t') 


-Our 


(108) 


/\  /S 


0 


This  expression  can  be  further  elaborated  on,  using  the  expression  of  the  ps-coupling 
in  terms  of  Coulomb  interactions.  It  tells  us  that  the  dipole  density  created  inside  the 
metal  at  distance  2,,  and  times  t’  is  propagated  in  the  s-region  and  perturbs  the  p-region 
at  times  t > t’  . 

In  the  dipole-dipole  approximation,  the  p-  and  s-  region  coupling  Hamiltonian  for 
the  electrons  and  nuclei  can  be  written  as 


0 


/S  /S 


Hps  — 


DpDsi^s ) 


dz 


z 


s 


(109) 


s 


— oc 


/N 

where  Dp  is  the  p-region  dipole  operator. 

Replacing  Eq.  (109)  into  expression  (108),  we  get  the  effective  dipole  expression  as. 


(f,  t']  2s)  = Dpf^{t',  t:  2s) 


(110) 


(Ill) 


— oc 


This  shows  that  within  the  first  order  perturbation  theory,  the  effective  primary  transition 
dipole  is  the  original  primary  dipole  modified  by  the  secondary  region  response,  but  it  is 
independent  of  the  electrical  field  strength  which  means  that  there  is  a linear  response  to 
the  field.  To  consider  the  nonlinear  response  case,  with  the  same  procedure,  we  include 

/N  ✓V 

r‘'(f)  to  all  the  orders  in  the  field  (or  fluence)  instead  of  just  keeping  the  first  term  Tq, 
and  we  can  write  the  effective  transition  dipole  as 


Ap[tj'\z_^,F)  = bpfs{t' 


(112) 


with  F the  fluence  value,  and 


0 


j\{t\t\Z!,,F)  =itr^{V-\t,F)[b^[zf,j')j,  j dz'. 


/ D,{z',J.)i 


.1 13 


]}/fr,[r'*(f)]  (113) 


's 


— oc 


It  is  obvious  that  the  transition  dipole  is  dependent  on  the  electric  field  strength  and  this 
gives  a nonlinear  response  to  the  field. 

As  an  extention  of  the  dipole-dipole  approximation,  the  p-  and  s-  region  coupling 
Hamiltonian  for  all  the  charges  can  be  generally  written  as 


/N 


Ppidp:  bPs{f}n:  f) 
Wp  ~ ^.s| 


(114) 


where  cjp  and  q,^  are  the  p-region  and  s-region  position  coordinates,  pp(qp,  t)  and  Ps{qs^  f ) 
are  the  p-region  and  s-region  charge  density  operators  respectively. 

Replacing  Eq.  (114)  into  effective  dipole  expression  (108),  we  get  the  effective 
dipole  expression  as. 


fbqppp{qp,t)f.,{f'j,: 


(115) 
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with 


J Wp  (/s 


Wp~<1p\"  ■ ■ (116) 

This  give  the  effective  transition  dipole  in  the  primary  region  for  a linear  response  to 
charge  density  fluctuations.  For  the  nonlinear  response,  similarly  to  the  discussion  in 
the  dipole-dipole  approximation,  we  can  write  the  primary  effective  transition  dipole  as 


Ap[tj':Zf.,F)  = J d^qppp{qpj,)f,^{t' ,t:  ,z,,,  F) 


(117) 


with 


\Qp  Qs 


(118) 


0 


A simplified  yet  meaningful  model  of  primary  effective  transition  dipole  can  be 
constructed  assuming  a fast  stochastic  response  for  both  linear  and  nonlinear  response 
cases,  giving 


Ap{t,  t';  z^,  F)  = 2S{t  - t')Ap{t:  Zg,  F) 


(119) 


At  an  average  depth  Zg  inside  the  metal;  the  perturbation  of  the  primary  region  would 


then  be  given  by 


0 


<?!,(<) 


dzg  I dt  2S{t  - t )Ap{t:  Zg,  F)Sioc{zg,  t) 


— oc  0 

0 


(120) 


= - / dzgAp{t]Zg,F)£ioc{zgJ) 


— OC 


and  is  written  in  terms  of  an  effective  primary  dipole  created  by  the  indirect  mechanism. 


CHAPTER  5 

NUMERICAL  APPROACH  FOR  WAVEPACKET  PROPAGATION 

Introduction 

Among  the  methods  for  quantum  mechanical  numerical  calculations,  the  ones  based 
on  time  dependent  wavepackets  have  now  become  important.  The  success  of  these 
methods  are  due  to  their  natural  correspondence  to  experiment  - starting  from  an 
initial  state,  following  the  event  through  time  and  generalizing  the  intuitive  classical 
description.  Moreover,  recent  progress  in  time  dependent  methods,  in  particular  the 
introduction  of  grid  representations  and  discrete  variable  representations,  have  made 
these  methods  competitive  with  other  exact  quantum  mechanical  calculations.  In  this 
chapter,  we  first  overview  general  time  dependent  propagation  methods,  then  present 
the  modified-split  operator  method  used  in  our  calculation  and  give  a brief  discussion 
of  how  to  combine  the  split-operator  method  with  a discrete  variable  representation  for 
problem  with  angular  degrees  of  freedom.  Finally  we  give  a brief  description  of  the 
computer  program  structure  and  its  implementation. 

Time  Dependent  Methods 


The  purpose  of  this  study  is  to  solve  the  time  dependent  Schroedinger  equation 


(121) 


with 


t 


F„  = Kj,  + U;,  - -Wj,  - Ap  • £{t) 


(122) 
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To  solve  this  equation  using  propagation  methods,  we  basically  need  two  steps:  the 
first  step  is  an  accurate  discrete  spatial  representation  of  the  wavefunction.  Once  the 
representation  is  constructed,  an  initial  wavefuction  can  be  propagated  in  time.  The 
second  step  is  to  choose  a propagation  methods. 

There  are  several  time  dependent  propagation  methods [62-64].  The  most  widely 
used  propagation  algorithms  in  quantum  dynamical  calculations  are:  the  Chebyshev 
polynomial  expansion  scheme [65-69];  the  second-order  differencing  scheme  [70];  the 
short  iterative  Lanczos  scheme  [63];  and  the  split  operator  propagation  scheme[49, 
71-74,  43].  These  algorithms  all  have  pros  and  cons,  and  are  used  for  various  classes 
of  dynamics  problem. 

/N 

In  this  study,  the  Hamiltonian  Fp  is  time  dependent  and  not  Hermitian  because 

/N 

it  includes  a dipole  coupling  term  • S{t)  and  a dissipation  term  Wp.  These  terms 
imply  that  we  can  not  use  the  first  mentioned  two  algorithms,  which  need  either  a time 
independent  Hamiltonian,  or  a strictly  Hermitian  Hamiltonian  for  an  efficient  and  stable 
scheme.  For  the  short  iterative  Lanczos  scheme,  the  evolution  matrix  depends  on  the 
initial  state,  which  makes  the  algorithm  complicated  when  we  need  to  calculate  the 
effects  of  different  initial  states. 

We  chose  the  split  operator  propagation  scheme  because  this  method  is  numeri- 
cally unconditionally  stable  and  it  is  very  efficient  for  problems  only  with  Cartesian 
coordinates.  For  every  time  step  the  computing  time  scales  as  NlogN,  where  N is  the 
number  of  grid  points.  After  some  modifications  of  the  original  split  operator  propaga- 
tion scheme,  it  can  treat  time  dependent  Hamiltonian  including  dissipation  in  Cartesian 
coordinates.  By  combining  the  modified  split  operator  algorithm  with  a discrete  variable 
representation  (DVR),  it  can  also  handle  rotational  degrees  of  freedom. 
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General  Initial  State 


The  initial  state  of  the  system  is  described  by  a two-component  wavefunction 


■0”(Q;  t=  0) 


'0g,v(Q;  0) 


0 


(123) 


This  initial  state  means  that  the  system  is  in  the  electronic  ground  state  before  the  light 
field  has  been  shining  on  the  adsorbate-substrate  complex.  However,  on  the  electronic 
ground  potential  surface  at  temperature  T,  the  initial  nuclear  vibrational  states  have  a 
Boltzmann  distribution  obtained  from 


(124) 


and  the  Boltzmann  weights  {/^.v}  are 


exp(-/).s6„^v) 

^exp(-/ii,,63,v) 


(125) 


where  /),  = 1/(^’bTs)  describes  the  dependence  of  the  initial  vibrational  state  distribu- 
tion on  the  temperature  T,  of  the  substrate. 


Modified  Split  Operator  Scheme  in  the  Electronic  Diabatic  Representation 


The  time  evolution  of  the  p-region  amplitude  for  nuclear  motions  can  be  obtained 
by  propagating  the  initial  wavepacket  in  time  increments  At,  as 


+ At))  = U^'^\t  + AtJ.)[iP^^\t)) 


(126) 


with  the  propagator  operator  given  for  small  At  by 


U^'^\t  + At,t)  = exp  [-iAtP^‘^\t)] 


(127) 


In  the  coordinate  Q representation,  it  gives 


(128) 
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Discretizing  it  in  space  grid  {Q„  } with  N points,  we  get 

IJ  At)  (129) 

with 

V^‘^\t+At,t)=  [(Q„/|(7('^>(t+Z\t^)|Q„,)]  (130) 

being  an  N x N matrix  and 

i’{t)  = [lP{Cln,t)]  (131) 

being  an  A/^  x 1 matrix.  We  can  also  discretize  the  time  variable,  with 

tn=tin+nAt  (132) 

If  we  assume  the  wavepacket  propagation  procedure  starting  from  An  = 0,  then  we  have 

tn=nAt  (133) 

After  discretizing  space  and  time  variables,  and  using  the  time-evolution  matrix  in  the 
diabatic  representation, 

Ai)=  exp[-?;z3fF^'^*(f„)]  (134) 

we  can  write  the  solution  formally  as 

fl— 1 

=n  U(‘^'(tn+i,A,)V’<'^^(f=0)  (135) 

n=0 

For  the  indirect  photodesorption,  we  have 

= K('^4v<'^4A(‘^)(f„)  • £{tn)  - ^W(^)(f„)  (136) 

If  we  define 

(f„)  - 


(137) 
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Eq.  (134)  can  be  approximated  by  the  split-operator  expression  accurate  to  the  second 
order[49] 


(138) 


xexp[—iAV^'^^{tn)/‘^]  + O(At^) 

Noticing  that  here  the  split  operator  scheme  is  extended  to  handle  time-dependent 
Hamiltonian  unlike  the  original  time-independent  Hamiltonian  scheme.  The  exponential 
of  the  complex  potential  energy  matrix  can  further  be  factored  into  exponential  for  its 
real  and  imaginary  parts, 


= exp[-AtW^'^^tn)/4]exp{-iAt[V^'^^+A^'’\tn)  ■ f(j!„,)]/2} 


xexp[—iAt^^^^ 


(139) 


xexp{-iAt[V^'^^+A^^^{tn)  • £{tn)]/2}exp[-AtW^^'{tn)/4]  + 0(A/") 

The  expression  (139)  for  the  time  evolution  over  time  dependent  coupled  complex 

energy  surfaces  can  be  called  the  “modified  split  operator  scheme”,  as  the  original 

split  operator  scheme  [75,  49]  is  employed  on  the  whole  propagator  to  separate  first 

the  imaginary  potential  contribution,  and  then  to  separate  its  kinetic  and  real  potential 

contributions.  We  note  that  both  and  • £{tn)  are  2 by  2 matrices 

with  off-diagonalize  elements,  and  we  need  to  diagonalize  them  before  we  can  do  the 

calculation.  If  we  define  the  real  potential  energy  part 


V'-(f„)  = V('^4A('^)(^„,).,f(f„)  (140) 


The  transformation  between  the  diagonalize  and  non-diagonalize  matrix  for  the  real 
potential  part  in  the  coordinate  representation  is  given  by 

Vh{tn)=  Lvv(t„)V'-(f.)L{,..(f„)  (141) 

and  the  transformation  matrix  Lyr  can  be  written  as 

cos((Q,f„)  .siii((Q,t„) 

-sinC(Q,t„)  co.sC(Q,^„) 


Lyr  (Q,  tf]  ) 


(142) 
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with  a rotation  angle  that  is  a function  of  the  p-region  coordinates  and  time, 


= -arctaii 


2V^(Q,t„) 


(143) 


Because  in  this  study  the  real  part  of  the  potential  is  time  dependent,  we  notice  that 
the  transformation  matrix  and  the  rotation  angle  will  also  be  time  dependent.  Due  to 
this  reason,  we  must  recalculate  them  for  each  propagation  time  step.  Similarly,  The 
imaginary  potential  is  diagonalized  into  W^)  by  a matrix  Lw  with  the  following 
relation, 

Wij(Q,f,)=Lvr-(QA)W('^>(Q,f,,)L{,,.(Q,t„,)  (144) 


where  we  also  assume  the  dissipation  is  time  dependent,  which  is  the  case  in  real 
photodesorption  processes.  The  transformation  matrix  Lpi’  can  be  written  as 


Lvv’(Q,  t,i) 


cos^(Q,t„)  sin^(Q,f„) 
-sin^(Q,f„)  cos^(Q,f„) 


(145) 


with  a rotation  angle  that  is  a function  of  the  p-region  coordinates  and  time. 


= - arctan 


2WciQ,tn) 


LW,(Q,f„)-W,,(Q,/,)J 


(146) 


Therefore  if  the  dissipation  is  time  dependent,  as  in  the  real  potential  part,  we  need  to 
recalculate  the  transformation  matrix  for  each  propagation  time  step.  With  these  trans- 
formations, the  propagator  matrix  of  the  time  dependent  Hamiltonian  with  dissipation 
can  be  written  in  the  form 


D 


X [l|;4  Q , )Uy.  (f„^  i , in  )Lyr  ( Q , /„ )] 


■D 


xU 


(d) 

K 


7),  )Lyr(Q,  f„)] 
77,)LrV’(Q,^77)] 


(147) 


51 


with 


j = exp{-iAtV^jj{tn)/2) 

(148) 

<’  = 

: exp 

(149) 

) = exp(-A^WD(*n)/4) 

(150) 

Finally,  the  kinetic  energy  matrix  does  not  satisfy  the  locality  condition  in  the 
coordinate  representation  because  of  its  second  derivatives  with  respect  to  coordinates 
Q.  However,  we  know  that  the  kinetic  energy  operator  is  local  in  the  momentum 
representation.  If  we  do  not  need  to  include  angular  coordinates  in  the  problem, 
then  we  can  use  the  fast  Fourier  transform  (FFT)  algorithm  [76]  to  transform  the 
wavepacket  from  coordinate  representation  to  moment  representation,  to  propagate  in 
moment  representation  and  then  to  transform  back  to  the  coordinate  representation.  The 
kinetic  energy  is  expressed  as  in  the  momentum  representation,  so  that 


(151) 


where 


Tk{Q)  = 


AQexp(27riA:?7),/A^)  0 

0 AQexp(27T7A;?7?,/A^) 


(152) 


is  a 2 X 2 matrix.  If  we  define 


= exp{-iAt 


(153) 


we  finally  obtain  the  propagator  matrix  with  all  the  matrix  elements  localized  and  all 
the  exponential  matrices  diagonalized. 
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(154) 


Angular  Degree  of  Freedom 


In  the  last  section,  the  introduced  split-operator  scheme  for  the  wavepacket  propa- 
gation works  efficiently  when  there  are  no  angular  coordinates.  But  for  problems  such 
as  desorption  coupled  to  a frustrated  rotation  vibrational  mode,  we  need  to  include  an- 
gular coordinates.  In  this  situation  we  can  not  use  the  fast  Fourier  transform  to  do 
transformation  between  moment  space  and  coordinate  space  which  is  needed  to  localize 
the  kinetic  energy  term.  Here  we  give  a brief  discussion  of  how  to  handle  this  problem 
by  combining  the  discrete  variable  representation[77-82]  into  the  split-operator  scheme. 

The  effective  Hamiltonian  of  the  p-region  is  given  by 


the  rest  which  only  depends  on  coordinate  degrees  of  freedom  and  time.  If  we  assume 


(155) 


and  dropping  the  p index,  we  can  rewrite  it  as 


F{tn)  — ^0  + P'li'f'-n) 


(156) 


where  Fq  is  the  part  of  the  total  Hamiltonian  with  known  eigenfunctions  and  Fi(f„)  is 


that  I'lpi),  with  i = 1,2, 3. ..77,  gives  the  eigenfunction  set  of  Fq,  then  we  have 


i'fpil'fpj)  — Fj 

N 


i=l 


(157) 
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where  Sij  is  the  Kronecker  delta  and  I is  the  identity  matrix.  The  above  expressions 
can  be  written  in  integral  form 


/r/Qi/;,:(Q)V>j(Q)=  6ij 

' N (158) 

/ r/Q^  0,:(Q)0,:(Q')  = 1 

?:=i 

Using  the  discrete  variable  representation  [77-79],  we  discretize  the  above  two  equations 


with  M quadrature  points. 


M 


^ Qrt  )WrtV’j(Qrt  ) = Sij 


(159) 


J]'0,:(Qa)(W«W,,)^'0,:(Qa)  = 


/=! 


where  are  the  quadrature  points  and  W„  the  corresponding  weights  with 


Qo  — ■ 'Qm)ci  and 


(fUi,  W2,  ■ ■ -W\f)a  respectively.  From  Eqs.  (159) 


we  can  define  an  77  x M transformation  matrix  T with 

Tia  = -0/(Qa)w| 


(160) 


where  N is  not  necessary  equal  to  M.  The  matrix  elements  of  Fq  and  in  the 

{|0j)}  basis  are  given  by 


(Fo)ij=  6ij{ipi\Fo[ipj)  = EiSij 


(161) 


and 


= f dQ'il;i{Q)F^{Qj.n)'ipj{Q) 

M 
a — 1 

It  turns  out  that 


(162) 


F,(^,0  =TlFf(f„)T 


(163) 
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with 


(164) 


an  M X M matrix. 


Combining  the  split  operator  scheme  with  the  discrete  variable  representation,  finally 
we  get  the  propagator  matrix 


(165) 


Notice  here  that  the  exponential  matrices  are  diagonalized  and  the  transformation  matrix 
can  be  easily  obtained. 


Based  on  the  time  evolution  operator  expression  (147)  constructed  from  a 
time— dependent  Hamiltonian  with  dissipation  and  two  coupled  complex  potential  energy 
surfaces  by  means  of  the  modified  split  operator  scheme  given  in  this  chapter,  we  have 
implemented  a Fortran  77  program  package  for  any  two-  dimensional  Cartesian  degree 
of  freedom  wavepacket  propagated  with  a time-dependent  Hamiltonian  and  a time- 
independent  dissipation.  This  program  is  an  extension  of  the  previous  time-independent 
Hamiltonian  wavepacket  propagation  program  for  CO  photodesorption  from  the  Ni(OOl) 
surface  by  Beksic  [83].  The  flow  charter  of  the  present  program  is  given  in  Appendix. 

The  main  program  (“photod.f ’)  consists  of  four  functional  parts  including  construct- 
ing the  initial  wavepacket  ( subroutine  “iwavepacket.f’),  constructing  potential  energy 
surfaces  and  their  couplings  (“pes.f ’),  propagating  the  wavepacket  one  time  step  (“prop- 
agate.f’)  and  calculating  properties  (“property. f’).  Evidently,  this  structure  assures  a 
great  flexibility  in  the  development  and  use  of  the  program.  For  example,  for  the  same 
phtodesorption  system,  different  initial  conditions  or  potential  energy  functions  can  be 


Computer  Program  Development 
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obtained  by  merely  modifying  the  corresponding  subroutines,  without  changing  the  rest 
part  of  the  program.  For  applications  to  other  photodesorption  systems,  we  only  need 
to  add  another  potential  energy  surface  subroutine  to  the  energy  surface  module  with- 
out making  any  other  changes  to  the  original  program.  This  makes  the  program  easily 
applicable  to  different  two  dimensional  photodesorption  problems.  In  this  study  we 
implemented  the  potential  energy  surfaces  and  their  couplings  for  CO  photodesorption 
by  femtosecond  light  pulses  from  Cu(OOl). 

Similarly  to  Beksic’s  [83]  program,  the  time  evolution  in  this  algorithm  is  sequen- 
tially done  in  time  steps  by  using  the  modified  split  operator  scheme  to  advance  by 
one  time  step  and  calculating  the  time  dependent  properties  every  few  time  step  until 
reaching  the  given  final  time.  The  time  advancement  is  executed  as  a sequence  of 
matrix  multiplications  given  by  Eq.  (147).  with  overwriting  of  the  wavepacket  after 
each  multiplication;  for  instance,  operating  with  the  real  potential  propagator  is  done 
as  a sequence 

x(nA)  = Lvv(nZ\)  • 
x(A,nA)  = Uyr(nA/2,  r?,A)  • 

(166) 

'0^'^^(A,nA)  = 

x(A,nA)  =Lyr(7?.A)  • ■i/;^'^*(A,  r?,A) 

■iJM\A,7iA)  = x(A,nA) 

Notice  we  have  to  recalculate  the  transformation  matrix  at  each  time  step  because  the 
Hamiltonian  is  time-dependent. 

The  output  of  the  time  evolution  program  consists  of  the  time  dependent  properties 
we  are  interested  in:  vibrational  distribution  in  the  two  electronic  states,  populations 
of  the  two  electronic  states  and  photodesorption  yields.  But  we  can  easily  add  other 
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properties  such  as  wavepacket  time  evolution  and  time  dependence  of  expectation  values 
of  the  primary  coordinates. 


CHAPTER  6 

APPLICATION  I:  INDIRECT  PHOTODESORPTION  OF  CO 
FROM  THE  Cu(OOl)  SURFACE  BY  FEMTOSECOND  PULSES 


Introduction 

With  the  theory  and  numerical  method  previously  developed,  we  consider  in  this 
chapter  the  indirect  mechanism  which  has  been  observed  experimentally  for  photons  off 
resonance  with  the  bond  excitation  energies  in  C0/Cu(001)[29,  39].  Our  model  assumes 
a CO  molecule  absorbed  with  its  axis  perpendicular  to  the  surface,  on  top  of  a Cu  atom. 
Our  reference  frame  has  its  origin  at  this  Cu  atom,  with  the  z-axis  perpendicular  to  the 
surface  toward  the  CO.  The  dynamics  of  desorption  is  assumed  to  involve  the  softest 
vibrational  mode  of  the  absorbed  CO,  which  is  its  frustrated  translation  mode,  or  “T 
mode”,  parallel  to  the  surface[21].  Therefore  one  coordinate  is  the  position  Z of  the 
center  of  mass  of  CO  above  the  surface  and  the  other  coordinate  is  taken  to  be  a 
displacement  x parallel  to  the  surface.  The  C — O distance  changes  from  its  value  when 
CO  is  bonded  to  the  surface  into  its  value  when  free  at  large  distances.  This  has  an 
effect  on  the  excitation  of  vibrations  of  CO  during  desorption.  Here  however  we  are 
focusing  on  the  desorption  rate,  which  depends  on  energy  transfer  between  translational 
and  vibrational  motions  of  the  center  of  mass  of  CO  and  should  not  be  very  sensitive 
to  the  CO  stretch  vibration  because  its  excitation  energy  is  much  larger  than  the  T- 
mode  excitation  energy  so  that  energy  is  transferred  mostly  between  the  T-mode  and 
translation.  We  take  the  C — O distance  to  be  constant  and  equal  to  its  value  dco  when 
CO  is  bonded  to  the  surface;  the  center  of  mass  of  CO  is  at  position  Z above  the  Cu 
atom,  and  the  C atom,  closest  to  Cu,  is  at  the  position  z.  The  displacement  of  the 
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CO  center  of  mass  parallel  to  the  surface  is  given  by  the  position  x.  We  therefore 
have  for  the  distance  from  C to  the  Cu  atom,  dccu  = and  the  relation 

z = Z — 'ydco  , 7 = MqI{Mc  + Mq).  The  coordinate  system  for  CO  vibrating  at 
the  Cu  surface  is  given  in  Fig.  (5). 

Construction  of  the  Potential  Energy  Surfaces  and  Their  Couplings 

For  the  real  part  of  the  potential  energy  surface,  we  can  perform  semiempirical 
molecular  orbital  calculations  using  parameterization  of  electron  integrals  with  incom- 
plete neglect  of  differential  overlap  appropriate  for  systems  containing  transition  metal 
ions  (the  ZINDO  parametrization)  [46-48],  For  the  CO/Cu(001)  case,  the  ground  and 
first  excited  electronic  states  of  CO/Cu(001)  have  been  modelled  by  clusters  of  CupCO, 
with  n = 1,6,  10.  The  ground  electronic  state  was  obtained  from  ZINDO/ 1 calculations, 
which  are  parametrized  to  provide  the  lowest  electronic  state.  It  is  known,  and  we  have 
found,  that  the  bonding  energy  and  equilibrium  distance  vary  with  the  size  of  the  cluster, 
although  they  remain  qualitatively  correct.  Since  their  experimental  values  are  known 
for  the  ground  electronic  state,  we  have  rescaled  our  curves  to  fit  experimental  values. 
Electronically  excited  states  were  obtained  from  ZINDO/S  calculations.  These  construct 
many-electron  states  as  linear  combinations  of  singly  excited  Slater  determinants  ob- 
tained from  a reference  determinant.  Excitation  energies  are  to  a good  approximation 
independent  of  cluster  size,  and  they  are  of  the  right  magnitude  even  when  a single 
Cu  atom  is  considered.  Taking  Cu6CO  as  an  example,  we  find  that  the  lowest  elec- 
tronic excitation  results  from  promoting  electrons  from  bonding  MOs  between  C and 
adjacent  Cu  atoms,  into  unoccupied  MOs  which  are  either  non-bonding  between  CO 
and  the  Cu6  cluster,  or  antibonding  between  them.  Electronic  rearrangement  changes 
the  magnitude  and  direction  of  the  dipole  of  the  whole  system,  indicating  that  there 
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is  electronic  density  transfer  between  CO  and  the  metal  atoms,  but  the  transfer  is  not 
well  described  by  the  frequently  assumed  picture  in  which  an  excited  metal  electron  is 
captured  into  a resonance  level  of  CO  to  form  a negative  ion.  In  a simplified  picture 
with  only  the  three  colinear  atoms  in  CuCO,  the  excitation  is  a combination  of  the 
types  7T  (7  and  tt  -7  8,  which  induces  a transition  dipole  perpendicular  to  the 
triatomic  axis,  and  therefore  parallel  to  the  surface.  For  the  ground  electronic  state,  the 
position  and  depth  of  the  attraction  well  were  scaled  to  agree  with  the  known  exper- 
imental values.  For  the  excited  electronic  state,  Saalfrank  and  Kosloff  [25]  point  out 
the  hot-electron  mediated  chemistry  seems  possible  only  for  not  too  large  excitation 
energies.  Further,  for  systems  with  several  photoactive  excited  states  it  seems  possible 
that  the  lower-lying  states  may  be  predominantly  populated  by  hot  electrons,  while  the 
higher  ones  are  accessible  only  through  direct  excitation.  Based  on  this,  the  excited 
electronic  state  potential  was  obtained  by  adding  the  calculated  excitation  energy  for 
the  lowest  excitation  with  a non-zero  transition  dipole,  to  the  ground  potential  energy. 
These  were  reproduced  by  potential  energy  surfaces  parametrized  with  Morse-like  po- 
tential functions  in  the  diabatic  representation  so  that  when  transformed  to  the  adiabatic 
representation  they  agreed  with  calculated  values. 

For  the  real  part  of  the  ground  PES  we  write 

Vf,{Z,x)  =Vj’{exp[-2o;3[[x2  {Z  - '^dcof]^  - Zg  -f  7f/co]] 

- 2 exp[-tt,[[.'r^  + {Z  - ^dcof]^  - Z^  -f  jdco]]}  (167) 

, McoiMZf  2 
2 

with 

i\{Z)  =:  - ha.7ih[2a^{Z  - Z^)]}  (168) 

where  Ilg(Z)  represents  the  Z-dependence  of  the  T-mode  frequency  of  CO  at  position  Z, 
is  the  T-mode  vibrational  frequency  of  CO  bonded  to  the  surface  at  the  equilibrium 
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point  Zg,  and  the  factor  between  brackets  decreases  this  frequency  down  to  zero  at  large 
Z;  the  distance  Zj  is  somewhere  larger  than  the  inflection  point  of  Vg.  Parameters  for 
the  ground  state  PES  given  by  Eq.  (167)  are  presented  in  table  (1) 

For  the  real  part  of  the  first  excited  PES  we  write 

Ve{Z,x)  =l/°(exp[-2a;4[.7;2  + (Z  - jdcof]'^  - + jdco]] 

- 2 exp  [-Op  [[•'5^^  + {Z  - Thof]^  - Ze+  'ydco]]}  + (169) 

, McoiUZ?  2 
2 

with 

ile{Z)  = ilj - ifa7?/?.[op(Z  - Zrf)]}  (170) 

where  represents  the  Z-dependence  of  the  T-mode  vibration  frequency  in  the 

electronically  excited  state,  (2e^  is  the  vibrational  frequency  of  CO  at  the  equilibrium 
point  Ze,  and  we  again  turn  off  the  frequency  at  large  Z.  Parameters  for  the  excited  state 
PES  given  by  Eq.  (169)  are  presented  in  table  (2). 

The  real  part  of  the  coupling  matrix  elements  of  the  potential  energy  matrices  in 
the  diabatic  representation,  labelled  by  an  index  “c”,  represents  the  diabatic  electronic 
coupling  between  the  two  electronic  states  and  can  be  significantly  different  from  zero 
only  near  the  crossing  of  PESs  at  short  distances;  it  is  written  as 

yO  j 

Vc{Z,x)  = -^{1  - tanh[ac[x^  + {Z  - '^dcof]^  - Zc  + jdco]}  (171) 

where  Z^  denotes  the  position  of  the  surface  crossing  with  x = 0.  Parameters  for  the 
real  coupling  term  given  by  Eq.  (171)  are  presented  in  table  (3) 

The  imaginary  parts  of  the  ground  and  excited  PESs  come  from  the  coupling  of 
the  CO  orbitals  to  the  intraband  and  interband  electronic  transitions  of  the  Cu  substrate. 


Table  1:  Parameters  for  the  CO/Cu(001)  electronically  diabatic  ground  PES. 
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Table  2:  Parameters  for  the  CO/Cu(001)  electronically  diabatic  excited  PES. 
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Table  3:  Parameters  for  the  CO/Cu(001)  diabatic  coupling  potential  energy. 
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respectively,  as  derived  in  Ch.  4.  These  couplings  must  vanish  for  large  values  of  Z, 
as  the  adsorbate  remains  asymptotically  in  its  ground  state.  For  them  we  use  the  forms 

W^{Z,  x)  = wjj  exp{-2«4[.7:'2  + {Z  - 'ydcof]^  - Zg  + 'jdco]}  (172) 

We(Z.x)  = exp{-2o!e[[.T^  + {Z  - 'ydcof]^  - Ze+  -ydco]}  (173) 

Here  the  dissipation  exponents  have  the  same  form  as  in  the  Morse  potentials. 
The  dissipation  strength  parameters  Wjj  and  Wg  originate  in  the  density  and  charge 
fluctuations  in  the  substrate  and  are  related  to  dissipation  decay  times. 

The  imaginary  part  of  the  coupling  matrix  elements  in  the  diabatic  representation 
represents  the  coupling  of  the  adsorbate  electronic  states  with  the  interband  electronic 
transitions  of  the  substrate,  and  like  the  imaginary  part  of  the  PES,  it  is  expected  to 
have  a short  range  contribution  coming  from  the  choice  of  electronic  representation. 
We  have  chosen  its  form  as 

Wc{Z,x)  = [W^{Z,x)  X Wg{Z,x)]^  (174) 

With  this  form,  our  choice  of  imaginary  potential  matrix  elements  gives  non-negative 
eigenvalues  for  dissipation  rates.  Parameters  for  the  imaginary  parts  and  their  coupling 
term  given  by  Eqs.  (172),  (173)  and  (174)  are  presented  in  table  (4).  Here  = h/xg 
corresponds  to  the  dissipation  width  of  in  the  ground  potential  surface  at  the  equilibrium 
position  Z = Zg  with  the  ground  electronic  state  life  time  t^,  and  = h.lTg 
corresponds  to  the  dissipation  width  of  the  excited  potential  surface  at  the  equilibrium 
position  Z = Zg  with  the  excited  electronic  state  life  time  Xg  much  shorter  than  that 
for  the  ground  electronic  state,  with  both  values  taken  from  experiments  as  shown  in 


what  follows. 


Table  4:  Parameters  for  the  CO/Cu(001)  imaginary  potential  terms. 
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The  transition  dipole  moment  between  the  adsorbed  electronic  state  of  CO  and  the 
desorbing  excited  state  was  also  calculated  and  is  a function  of  the  adsorbate  nuclear 
coordinate  only,  given  by 

D{Z,x)  = DqZ^{1  — ta.7ih[T]{Z  — Z^)]}  (175) 

which  fits  its  calculated  values.  Its  parameters  are  presented  in  table  (5). 

To  minimize  the  number  of  parameters  in  these  functions,  we  have  used  the 
relations  , ac  = , and  at  the  crossing  point  of  the 

two  PESs.  Although  these  choices  are  somewhat  arbitrary,  they  give  a correct  order 
of  magnitude;  the  calculated  results  are  not  sensitive  to  the  accuracy  of  these  values 
because  the  potentials  crossing  is  located  at  a small  Z,  where  wavepackets  have  very 
small  amplitudes. 

In  addition  to  the  parameters  fit  to  calculations,  the  experimental  values  used  in  the 
PESs  are  Zg  = 3.59  au  and  rco  = 2.17  au  [84],  Vg®  = 0.0268  au  = 0.73  eV  [85],  Og 
= 1.52  au  coming  from  u;t(Cu-CO)  = 345  cm~'  [86],  and  fl^g  =1.45x10“^  au  = 3.94 
meV[87].  For  the  strength  of  the  ground  dissipative  potential  we  used  the  experimental 
rate  (inverse  of  the  relaxation  time)  known  from  experiments  on  vibrational  relaxation,  so 
that  = hlr^  comes  from  Xg  = 5 ps[33].  For  the  strength  of  the  excited  dissipative 
potential  we  used  the  relation  = hlxf  where  the  lifetime  of  the  electronically 
excited  states  is  believed  to  be  in  the  range  10“^  to  10*  fs  [88-91],  we  chose  Tg  =5, 
25,  50  fs.  These  parameter  values  are  close  to  the  expectation  values  of  the  dissipative 
functions  averaged  over  the  nuclear  wavefunctions,  and  describe  the  dissipation  effects 
that  electronic  excitations  in  the  substrate  have  on  vibrational  motions  during  desorption. 


Table  5:  Parameters  for  the  CO/Cu(001)  transition  dipole  function. 
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The  real  PESs  (167)  and  (171),  as  functions  of  the  adsorbate-substrate  distance  Z, 
for  a fixed  value  of  x = 0 corresponding  to  the  ground  electronic  state  equilibrium 
frustrate  vibrational  state  are  illustrated  in  Figure  (6). 

Form  of  the  Light  Pulses 

Although  the  theory  and  algorithm  developed  in  this  study  are  applicable  for  any 
form  of  light  pulse,  the  electric  field  pulse  is  chosen  to  have  a Gaussian  envelope.  To 
make  our  work  comparable  with  the  known  experimental  results,  the  carrying  wave  has 
a central  wavelength  Aq=620  nm  and 

S{t)  = ]co.s(o;o  • t)  (176) 

where  ujq  = ^ is  the  central  frequency  of  the  pulse  describing  the  rapid  oscillation 
of  the  light  field,  A is  the  temporal  width  of  the  pulse  [full  width  at  half  maximum 
(FWHM)]  and  F its  fluence  averaged  over  frequency  oscillations.  We  have  chosen 
A = f,Q.  With  this  pulse  form,  we  will  give  a detailed  discussion  of  photodesorption 
yields  of  CO  from  metal  surface  Cu(OOl)  dependent  on  properties  such  as  pulse  fluence, 
width  and  frequency. 

Numerical  Parameters  for  Wavepacket  Propagation 

There  are  two  critical  parameters  for  numerical  wavepacket  propagation,  the  time 
step  and  grid  spacings.  With  small  time  step  and  grid  spacings,  we  can  usually  get  a 
more  accurate  numerical  results  but  it  will  also  need  more  computer  time  and  memory 
space. 

For  the  grid,  special  requirements  appear  because  a fast  Fourier  transformation  algo- 
rithm was  included  in  our  split  operator  scheme  for  transformation  between  configuration 
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Figure  6:  Diabatic  potentials  for  ground  and  excited  electronic  states  versus  the  distance  Z from  the 
Cu  atom  at  the  surface  and  origin  of  coordinates,  to  the  center  of  mass  of  the  CO 
molecule  located  on  top  and  with  its  axis  perpendicular  to  the  surface.  The  CO 
center  of  mass  is  at  its  equilibrium  position  x=0  along  an  axis  parallel  to  the  surface. 
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space  and  momentum  space  representations.  For  an  accurate  representation  of  the  kinetic 
energy  K = ^ maximal  allowable  grid  spacings  are  estimated  by 

AZ  <<  2'KjPz 


A,t  <<  2tt  / Px 


(177) 


which  implies  choosing  about  15  grid  points  per  unit  wavelength. 

For  the  time  step,  and  although  the  split  operator  method  is  unconditionally  stable 
and  it  allows  one  to  choose  the  time  step  according  to  the  energy  range  of  interest,  a 
small  time  step  is  still  needed  because  the  Hamiltonian  is  time-dependent  with  a highly 
oscillating  electronic  dipole  coupling  term.  After  many  test  calculation,  for  a pulse  with 
wavelength  around  620  nm  and  period  of  2 fs,  we  give  the  numerical  chosen  parameters 
of  wave  propagation  in  table  (6). 

Time-dependent  Properties 

For  the  propagation  of  the  wavepacket,  starting  with  the  primary  region  in  its  ground 
vibrational  state  and  since  the  Z and  x variables  are  normal  coordinates  in  the  ground 
electronic  PES,  we  have  for  the  initial  packet 

\^UZ,x,0))  = ^,o{Z,x)\rJ{Z,x)) 

^go{Z,x)  = Aexp[— G!|-(Z  — exp[— a^.7;^]  (178) 

a|'  = 40.31  au  , a].  = 3.78  au 

To  complete  the  specification  of  input  data,  we  need  the  effective  primary  dipole  which 
follows  from  our  treatment  of  the  indirect  excitation  process.  This  is  the  most  difficult 
quantity  to  specify,  in  the  absence  of  a model  of  the  electronic  coupling  between  p- 
and  s-regions.  We  expect  it  to  be  proportional  to  the  primary  dipole.  In  effect,  the 
effective  dipole  involves  the  Coulomb  coupling  Hpg  of  the  p-  and  s-regions  as  found  in 
Eq.  (108);  an  expansion  of  this  factor  in  multipoles  would  show  a leading  term  coming 


Table  6:  Grid  parameters  for  the  two-dimensional  wavepacket  propagation. 
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from  the  interaction  of  the  primary  and  secondary  dipoles,  with  the  latter  averaged  with 
the  secondary  density  operator.  In  the  calculations  that  follow  we  therefore  make  a 
choice  with  the  effective  transition  dipole  in  the  primary  region  written  as 

Lp  = bp{a  + (179) 

where  Dp  is  the  transition  dipole  from  the  ZINDO  calculations,  a and  fi  are  parameters 
which  can  be  obtained  from  experiments,  F is  the  pulse  fluence,  and  Fq  = 1.0  mj/av?. 
With  n equal  to  several  values,  we  can  calculate  the  effect  of  nonlinear  transition  dipoles 
on  photodesorption. 

Results  for  Non-linear  Response 

Figure  (7)  and  (8)  present  the  comparison  between  one  and  two  dimensional 
calculation  results  for  the  pulse  fluences  of  1.5  mJ/cm^  and  3 mJ/cm^  respectively. 
In  both  case  the  one  dimensional  calculations  give  about  7%  more  yields  than  the  two 
dimensional  calculations.  These  can  be  explained  in  that  part  of  the  energy  goes  to  the 
excitation  of  CO  frustrated  translational  vibrations.  The  two  dimensional  calculations 
are  more  realistic  than  the  one  dimensional  calculations. 

Figure  (9)  through  (12)  present  the  time  evolution  of  electronic  ground  and  excited 
states  populations  for  both  one-  and  two-  dimensional  calculations.  There  are  rapidly 
oscillations  during  pulse  duration  which  means  that  there  exist  multiple  excitation  and 
deexcitation  steps.  The  difference  between  the  excited  asymptotic  populations  for  one 
and  two  dimensional  calculations  are  around  1%. 

Figure  (13)  present  desorption  yields  dependence  on  pulse  fluence  for  both  one 
and  two  dimensional  calculations.  In  both  calculations  the  yields  show  nonlinear 
dependence  on  the  pulse  fluence.  From  the  comparison  between  the  one  and  two 
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Figure  7:  Time  evolution  of  photodesorption  yields;  comparison  between 
one  and  two  dimensional  calculations  for  pulse  fluence  of  1.5  mJ/cm^ 
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one  and  two  dimensional  calculations  for  pulse  fluence  of  3.0  mJ/cm 
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between  one  and  two  dimensional  calculations  for  pulse  fluence  of  1.5  mJ/cm 
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between  one  and  two  dimensional  calculations  for  pulse  fluence  of  1.5  mJ/cm 
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between  one  and  two  dimensional  calculations  for  pulse  fluence  of  3.0  mJ/cm 
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between  one  and  two  dimensional  calculations  for  pulse  fluence  of  3.0  mJ/cm 
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dimensional  calculations,  we  have  seen  that  both  give  similar  desorption  results  within 
a 7%  difference. 

Figure  (14)  and  (15)  present  population  dependences  on  pulse  fluence  for  both 
one  and  two  dimensional  calculations.  In  both  calculations  the  populations  also  show 
nonlinear  dependence  on  the  pulse  fluence.  From  the  comparison  between  the  one  and 
two  dimensional  calculations,  we  have  seen  that  both  give  similar  population  results 
within  a 7%  difference.  From  now  on  in  this  chapter  all  the  results  are  from  one 
dimensional  calculations. 

Figure  (16)  presents  the  time  evolution  of  photodesorption  yields  with  several  pulse 
fluences.  It  shows  that  desorption  happens  very  quickly  with  time.  Figure  (17)  and  (18) 
present  the  electronic  ground  and  excited  states  population  time  evolution  respectively. 

Figure  (19)  presents  the  photodesorption  yield  dependence  on  pulse  widths.  Figure 
(20)  presents  the  photodesorption  yield  dependence  on  pulse  carrying  wavelengths. 
Figure  (21)  presents  the  photodesorption  yield  dependence  on  substrate  dissipation 
time  constants.  Theses  figures  show  that  the  yields  are  strongly  dependent  on  carrying 
wavelength  and  dissipation  rates,  but  not  much  on  pulse  widths. 

Figure  (22)  through  Figure  (27)  present  the  excited  and  ground  amplitude  distribu- 
tion along  coordinate  Z at  different  times.  The  total  amplitude  starts  as  a simple  Gaussian 
function  on  the  ground  electronic  state  at  the  initial  time  t=0;  as  time  increases,  the  am- 
plitude maximum  moves  from  its  initial  positions,  and  also  simultaneously  change  its 
shape.  Finally  the  amplitude  splits  into  two  parts,  with  a very  small  part  desorbed  and 
the  major  part  remaining  adsorbed  at  the  surface. 
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Figure  13:  Photodesorption  yields  dependence  on  pulse  fluence,  for  one  and  two  dimensional  calculations 
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fluence;  comparison  between  one  and  two  dimensional  calculations 
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fluence;  comparison  between  one  and  two  dimensional  calculations 
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Figure  19:  Time  evolution  of  photodesorption  yield  dependence  on  pulse  wi 
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Figure  20:  Time  evolution  of  photodesorption  yield  dependence  on  pulse  central  wavelengths 
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Figure  21:  Time  evolution  of  photodesorption  yield  dependence  on  substrate  dissipation 
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Figure  25:  Ground  amplitude  distribution  along  the  coordinate  Z at  the  times  t=0  au,  t=250  au  and  t=500 
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ure  26:  Ground  amplitude  distribution  along  the  coordinate  Z at  the  times  t=1000  au, 
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Results  for  Chirping 

We  now  consider  the  pulse  chirping  effects  on  photodesorption.  The  pulse  is  still 
chosen  to  have  a Gaussian  envelope,  and  a carrying  wave  with  a central  wavelength 
Ao=620  nm,  but  now  the  frequency  is  linearly  dependent  on  time  with  the  following 
form, 


£{t)  = 2F^/2(-^)^/'‘exp[-^^^-^^^]cos[a;(f)  -t] 
Lu{t)  = ojQ  at 


(180) 


with  ujQ  = « = p'^,  where  p is  a chirping  parameter,  A the  pulse  width  , and  F 

its  fluence  averaged  over  frequency  oscillations.  We  have  chosen  A = to  = 100  fs  as 
before.  For  the  effective  transition  dipole,  we  use  the  same  form  and  parameters  as  in 
the  previous  calculations  in  this  chapter. 

Figure  (28)  and  Figure  (29)  present  the  negative  and  positive  chirping  effect  on  the 
electronic  ground  state  population,  respectively.  Note  that  negative  and  positive  chirping 
have  different  effects  on  the  population,  comparing  to  zero  chirping  pulse:  the  negative 
chirping  decreases  the  ground  population  and  the  positive  chirping  increases  the  ground 
population.  Figure  (30)  and  Figure  (31)  present  the  negative  and  positive  chirping  effects 
on  the  electronic  excited  populations,  respectively.  As  before,  they  show  different  effects 
on  excited  populations:  the  negative  chirping  increases  the  excited  populations  and  the 
positive  chirping  decreases  the  excited  populations.  Notice  that  the  sum  of  ground  and 
excited  population  is  not  unity  because  of  the  dissipation  effects  of  the  substrate.  Figure 
(32)  shows  the  negative  and  positive  chirping  effects  on  photodesorption  yield.  From 
the  figure,  we  can  see  that  the  negative  chirping  increases  the  photodesorption  yields 
and  the  positive  chirping  decreases  the  photodesorption  yields,  respectively.  This  effect 
is  quite  pronounced  and  indicates  the  possibilities  of  optical  control. 
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Notice  also  that  all  the  populations  oscillate  rapidly  before  finally  coming  to  an 
asymptotic  value  which  depends  on  the  transition  dipole,  pulse  fluence  and  substrate 
dissipation.  But  the  photodesorption  yields  keep  increasing  smoothly  with  time  until 
finally  reach  an  asymptotic  value  which  depends,  of  course,  on  those  quantities. 

Comparison  with  Experiment  Results 

We  now  compare  our  calculations  with  the  experimental  results  presented  by  Prybyla 
et  al.  [29].  The  experiments  were  performed  on  CO/Cu(lll)  system  at  coverage  of 
0.52  monolayer(ML)  at  a temperature  of  95  K with  100  fs  laser  pulses  operating  at 
around  620  nm.  They  found  that  CO  molecules  desorb  from  a Cu(lll)  surface  within 
325  fs  after  irradiation  with  an  100  fs,  2 eV  photon  energy  pulses.  They  also  found 
that  the  desorption  yield,  which  was  measured  with  a mass  spectrometer,  increased 
nonlinearly  with  the  intensity  and  can  be  fitted  by  the  power  law  Y = const  For 
a pulse  fluence  of  4.5  mJ/cm^,  the  yield  is  0.05  ML.  In  our  calculation,  we  choose 
a Gaussian  pulse  with  the  same  properties  as  the  experiment.  We  further  choose  the 
effective  transition  dipole  Ap  as 

Ap  ^ Dp{a  + (181) 

where  Dp  is  the  calculated  primary  transition  dipole  from  the  ZINDO  calculation,  and 
Fq  = 1.0  mj/cnr.  In  this  model,  we  let  n equal  2.  Then  by  fitting  the  calculated 
yields  for  F = 1.5  inJ /cmr  and  F = 3.0  nij/cm?  to  experiment,  we  can  get  a — 0.4 
and  P = 0.28.  With  these  parameters,  figure  (33)  shows  the  calculation  yields  and 
experiment  yield  from  Prybyla  et  al..  For  all  the  measured  fluences,  from  the  figure  we 
can  see  that  our  results  are  very  consistent  with  the  experimental  results  for  fluences  in 
the  range  of  F = 1.0  7iiJ/cnP  to  F = 5.0  mJ/aiP. 
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Figure  28:  Negative  chirping  effect  on  the  population  of  the  ground  potential  surface 
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Figure  31:  Positive  chirping  effect  on  the  population  of  the  excited  potential  surface 
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Figure  32:  Positive  and  negative  chirping  effects  on  photodesorption  yields 
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Figure  33:  Photodesorption  yield  dependence  on  fluence;  comparison  with  experiment 
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Now  we  move  on  to  calculate  the  desorption  time  with  the  same  parameters.  Figure 
(34)  shows  a Gaussan  pulse  centered  at  100  fs,  and  the  corresponding  photodesorption 
yield  time  evolution.  From  the  figure,  we  can  see  that  the  pulse  starts  approximately  at 
t = 0.0  fs  and  the  photodesorption  process  is  completed  around  t — 280  fs,  and  this  is 
quite  consistent  with  the  experiment  result  from  Prybyla  et  al.,  that  t is  less  than  325  fs. 
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Figure  34:  Time  evolution  of  photodesorption 


CHAPTER  7 

APPLICATION  II:  ISOTOPE  AND  TEMPERATURE  EFFECTS  ON  THE 
INDIRECT  PHOTODESORPTION  OF  CO  FROM  THE  Cu(OOl)  SURFACE 


Introduction 

As  a second  application  of  the  theory  developed  in  this  thesis,  this  chapter  presents 
the  results  on  the  dependence  of  the  desorption  yields  on  the  CO  isotope  masses  and  on 
the  initial  substrate  temperature.  A brief  discussion  is  given  to  compare  the  calculations 
with  the  known  experimental  results. 

Photodesorption  Yields  Dependence  on  Isotopes 

In  this  section  we  study  the  dependence  of  CO  photodesorption  yields  on  the 
isotopic  substitution  of  the  CO  adsorbate.  We  know  that  within  the  Born-Oppenheimer 
approximation,  different  isotopes  of  a given  atom/molecule  move  on  the  same  potential 
surfaces.  But,  because  of  changes  in  the  nuclear  masses,  the  dynamics  of  the  system 
changes  with  isotopic  substitution  and  causes  photodesorption  yield  changes.  A formula 
was  derived  with  the  assumption  that  the  primary  excitation,  the  potential  surface  in  the 
excited  state,  and  the  lifetime  of  the  excited  state  are  independent  of  the  isotopic  mass, 
and  gave  the  dependence  of  the  photodesorption  yield  on  the  mass  of  the  particle[92], 

Y = Aexp(—cM2)  (182) 

where  A and  c are  constant  for  one  adsorbate/substrate  system  and  M is  the  mass 
of  the  desorbing  particle.  That  assumption  therefore  gives  an  isotope  effect  for  the 
photodesorption  yield:  the  light  isotopes  escape  more  rapidly  from  the  surface  and. 
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consequently,  are  expected  to  give  large  photodesorption  yields  than  the  heavy  ones. 
This  isotope  effect  was  observed  in  several  system  [93,  94]. 

Recently,  Szabo  and  Yates  [95]  reported  an  anomalous  desorption  yield  in  an  elec- 
tron stimulated  desorption  experiment  of  CO/Pt(l  1 1)  system.  They  found  that  instead  of 
the  light  isotopes  having  higher  yields,  the  desorption  yields  always  followed  the  trend: 
yC^C^^O)  > yC^C^^O)  ~ ^ y(i3ci6o)  independent  of  the  CO  coverage 

and  primary  excitation  energies.  More  recently,  for  the  photodesorption  of  CO  from 
Cu(l  1 1)  system,  Tom  et  al.  [96]  found  that  the  yields  for  the  heavier  masses  are  always 
higher  than  for  the  lighter  masses  for  the  3 isotopes,  ^^C^^O,  ^^C^^O  and  ^^C^*^0,  they 

10  1 A 

tried.  They  did  not  do  C O because  this  mass  is  experimentally  ambiguous  with 
N2.  They  found  that  the  isotopic  differences  vary  slightly  with  coverage  but  the  trend 
is  always  there  and  gives  Y(^^C^®0)  > y(^^C^®0)  > y(^^C^^O). 

Here  we  present  our  results  with  four  combinations  of  C and  O isotopes  using  the 
model  we  developed  in  this  thesis  with  the  isotope  parameters  present  in  table  (7),  in 
atomic  units. 

Figure  (35)  shows  a two  dimensional  calculations  of  the  photodesorption  yields 
dependence  on  isotopes.  The  parameters  for  these  calculation  were  chosen  the  same 
as  in  the  tables  of  chapter  6 with  the  isotopes  parameters  given  in  table  (7).  The 
calculations  clearly  show  an  anomalous  yields  dependence  on  isotope  masses,  with 
y(^^C^^O)  > which  is  in  good  agreement 

with  the  experimental  trend  given  by  [95].  But  our  results  are  different  from  those 
from  Tom  et  al.[96]. 

To  find  whether  the  anomalous  desorption  yields  are  due  to  the  frustrated  translation 
motion  or  to  the  center  of  mass  change  of  the  different  isotope  molecules,  we  did 
one-dimensional  calculations  with  the  same  parameters  as  in  the  two-dimensional 


Table  7;  Parameters  for  the  CO  isotopes  (masses  and  distanees  in  atomic  units) 
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Figure  35:  Two-dimensional  calculations  of  the  photodesorption  yields  dependence  on  isotope  masses 


no 


calculations,  except  that  we  restricted  the  CO  motion  to  one  dimension.  Figure  (36) 
shows  the  one-dimensional  calculations  of  the  photodesorption  yields  dependence  on 
isotopes.  It  is  clear  that  the  anomalous  desorption  yields  still  hold  as  in  the  two 
dimensional  calculations  but  with  smaller  yield  differences  for  each  isotope.  It  seems  this 
anomalous  desorption  yield  phenomenon  comes  mainly  from  center  of  mass  differences 
of  the  CO  isotopes.  From  table  (7),  it  is  obvious  that  the  desorption  yields  depend  on 
the  value  of  7,  with  712  > 713  q > 712  (^ir,  q > 71,3  q.  This  is  also  in  good 

agreement  with  the  qualitative  explanation  given  by  Szabo  and  Yates  [95]. 

Figures  (37)  and  (38)  show  the  effects  of  changing  the  pulse  carring  frequency 
on  the  photodesorption  yields  for  several  isotopes.  Increasing  or  decreasing  the  pulse 
frequency  produces  some  yield  changes  for  each  isotope,  but  the  anomalous  desorption 
yield  dependence  still  stays  the  same. 

Figures  (39)  and  (40)  show  the  substrate  dissipation  effects  on  the  desorption  yields 
for  several  CO  isotopes.  With  the  lifetime  for  the  excited  potential  surface  change  as 
Te=25  fs  and  Tg=5  fs,  the  results  give  the  same  desorption  trend  as  the  calculation  with 
Te=50  fs,  but  with  a much  smaller  desorption  yield  for  each  CO  isotope  when  Tq=5  fs. 

Figures  (41)  and  (42)  show  the  potential  surface  effects  on  the  photodesorption 
yields  for  several  isotopes.  In  this  calculation,  we  shift  the  excited  potential  surface 
either  inward  (toward  the  surface)  with  Ze=3.90,  or  outward  with  Ze=3.94.  The  results 
show  that  the  desorption  yields  are  very  sensitive  to  the  potential  surface  changes,  but 
the  same  desorption  trends  still  hold  for  the  isotope  effects. 

Photodesorption  Yields  Dependence  on  Initial  Substrate  Temperature 

In  this  section  we  study  the  photodesorption  yields  dependence  on  the  initial 
substrate  temperature.  In  most  experiments  the  substrate  temperatures  are  around  100  K. 
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Figure  36:  One-dimensional  calculations  of  the  photodesorption  yields  dependence  on  isotope  masses 
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for  several  isotopes:  calculations  with  wavelength  A=600  nm 
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for  several  isotopes:  calculations  with  wavelength  A=640  nm 
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This  temperature  is  not  enough  to  excite  the  C-Cu  bond,  but  it  is  large  enough  to  excite 
the  frustrated-translation  vibrations.  Similar  to  Eq.  (178),  the  initial  packet  including 
vibrational  excitation  in  x-direction  is  given  by 


n 


^.70vx(2’,  x)  = A^exp[— o;|(2'  - Zg)^]  exp[-alx^]H^^{'/2a3.x)  (183) 


2 


= 3.78  av, 


o;^  = 40.31  fl/u  , 

where  Hy^{\/2axx)  is  the  -w^rth  order  Hermite  polynomial  whose  values  for  from 


0 to  2 are 

Ho{V2oixx)  = 1 

Hi{\/2axx)  = 2\/2axX  (184) 

H2{\/2axx)  = —(1  — 4a^x^) 

Here  we  give  the  results  for  the  T-vibrational  motion  initially  in  the  first  and  second 
excited  states.  In  Figures  (43)  and  (44),  state-resolved  yields  Yv^{t)  are  given  for 
Vx  = 0,  1 and  2.  The  results  show  how  vibrational  excitation  causes  the  desorption 
yields  to  increase,  because  excited  wavepackets  occupy  a larger  fraction  of  phase  space, 
and  have  wider  energy  distributions.  Therefore,  our  frustrate-translation  model  predicts 
that  surface  heating  should  increase  the  desorption  yields  since  the  T-vibrations  can  be 
easily  excited. 

In  summary,  we  have  seen  in  this  chapter  that  our  calculations  are  in  agreement  with 
the  known  experimental  results  [95,  96]  for  the  CO  isotope  effects  on  desorption  yields, 
which  show  an  anomalous  desorption  yields  dependence  on  the  CO  isotope  masses.  It 
seems  that  the  anomalous  yields  are  coming  from  changes  of  distance  of  their  center 
of  mass  to  the  Cu  surface.  We  have  also  seen  that  within  the  frustrate  translational 
desorption  model,  the  T-vibrations  can  be  easily  excited  and  therefore  the  desorption 
yields  are  sensitive  to  the  initial  substrate  temperature. 
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Figure  43:  Temperature  effects  on  the  photodesorption  yields: 
the  initial  T-vibration  in  the  first  excited  and  ground  states 
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CHAPTER  8 

CONCLUSIONS  AND  DISCUSSION 


Contributions  to  the  Theory  of  Femtosecond  Photodesorption 

A theory  has  been  developed  for  direct  and  indirect  femtosecond  photodesorption. 
The  theoretical  method  developed  in  this  dissertation  describes  the  desorption  induced  by 
multiple  electronic  transition  (DIMET),  caused  by  an  interaction  of  a femtosecond  UV 
light  pulse  with  the  adsorbate-substrate  complex.  Our  treatment  has  been  based  on  the 
partition  of  the  whole  system  into  primary  and  secondary  regions,  and  the  derivation  of  a 
quantum  Langevin  equation  assuming  a stochastic  secondary  region  which  may  interact 
with  both  primary  region  and  a femtosecond  light  pulse.  In  the  present  work  we  have 
elaborated  on  the  indirect  photodesorption  process,  where  a local  electromagnetic  field 
in  the  substrate  induces  an  effective  primary  dipole. 

We  derived  equations  for  this  effective  primary  dipole  with  a treatment  of  the 
response  of  the  secondary  region  to  the  local  field.  The  field  creates  electronic  excitations 
in  the  substrate  which  couple  through  Coulomb  interactions  to  the  primary  region,  to 
produce  an  effective  transition  dipole  there.  Within  perturbation  theory,  we  found  that 
to  first  order  the  primary  transition  dipole  will  be  linearly  enhanced  by  the  light  pulse; 
to  higher  order,  we  found  it  will  be  nonlinearly  dependent  on  the  light  pulse  fluence. 

The  energy  dissipation  term  appears  in  the  primary  region  Hamiltonian  as  an 
anti-Hermitian  addition  to  the  potential  energy  function.  Its  form  was  derived  for 
electron-hole  pair  excitations  with  the  knowledge  that  lattice  vibrational  excitation  is 
a much  slower  process.  The  form  of  the  dissipative  potential  depends  upon  a time 
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correlation  function  of  the  charge  density.  The  self-consistent  field  coupling  of  the 
primary  and  secondary  region  is  crucial  to  the  derivation  of  the  dissipative  potential, 
as  the  charge  density  fluctuation  time  correlation  function  vanishes  in  its  absence.  The 
energy  dissipation  in  desorbing  species  dynamics  is  a direct  consequence  of  the  self- 
consistent  field  treatment  of  the  system. 

The  experimental  conditions  of  photodesorption  with  electronic  transitions  may 
involve  either  the  adsorbate  interacting  with  a short  pulse  of  light  directly  and  inducing 
a transition  to  the  excitation  state,  or  the  light  pulse  can  first  interact  with  the  substrate 
creating  electron-hole  excitations  and  then  transfer  energy  to  the  primary  region.  By 
introducing  an  effective  transition  dipole  in  the  primary  region,  we  can  handle  both 
direct  and  indirect  photodesorption  by  using  the  same  propagation  procedure. 

Computational  Aspect  of  CO  Photodesorption  from  the  Cu(OOl)  Surface 

The  CO  femtosecond  photodesorption  dynamics  was  studied  by  solving  the  time 
dependent  Schroedinger  equation  with  the  time  dependent  effective  Hamiltonian  of  the 
primary  region  of  the  CO/Cu(001)  system.  The  desorption  was  modelled  as  a process 
with  two  electronic  states,  and  the  two  electronic  states  were  described  by  two  time 
dependent  complex  potential  energy  surfaces  so  the  effects  of  the  energy  dissipation 
in  the  substrate  and  the  interaction  between  light  pulse  and  substrate  were  included 
from  the  outset  electronic  structure  calculations,  based  on  a semi-empirical  molecular 
orbital  approximation  suitable  for  compounds  of  transition  metals,  were  carried  out  on 
Cu  clusters  to  provide  information  on  the  ground  and  excited  potentials  and  transition 
dipoles  for  CO/Cu(001).  An  analysis  of  the  ground  and  excited  electronic  states  shows 
that  there  is  extensive  electronic  rearrangement  during  excitation  due  to  transfer  of 
electrons  from  bonding  to  non-bonding  and  antibonding  orbitals  in  the  primary  region. 


123 


The  charge  redistribution  does  not  however  indicate  that  an  electron  is  captured  by 
the  absorbed  CO,  as  sometimes  postulated.  The  potential  energy  surfaces  and  their 
couplings  were  parametrized  semiempirically,  using  the  available  experimental  data  for 
the  CO/Cu(001)  system  and  electronic  structure  calculations  for  CO  bound  to  a single 
Cu  atom  and  an  Cu^  cluster. 

Our  model  dealt  with  a CO  molecule  absorbed  with  its  axis  perpendicular  to  the 
surface,  on  top  of  a Cu  atom.  The  dynamics  of  desorption  was  assumed  to  involve  the 
softest  vibrational  mode  of  the  absorbed  CO,  which  is  its  frustrated  translation  mode, 
or  “T  mode”,  parallel  to  the  surface.  Potential  energy  functions  and  couplings  were 
parametrized  in  a diabatic  representation,  using  calculated  and  measured  properties  of 
the  adsorbate. 

The  primary  time  dependent  Schrodinger  equation  was  solved  by  propagation 
of  wavepackets,  using  the  split-operator  propagation  scheme  to  construct  the  time 
evolution  operator,  with  a modification  to  include  the  imaginary  dissipative  potentials 
and  the  time-dependent  interaction  with  light.  The  diabatic  electronic  representation 
was  used  in  the  time  propagation,  because  it  eliminates  momentum  couplings  between 
the  electronic  basis  states;  the  real  and  imaginary  potential  energy  matrices  were 
diagonalized  separately  by  orthogonal  transformations  in  each  time  step.  The  time 
propagation  algorithm  was  numerically  implemented  using  a discrete  grid  to  represent 
the  coordinate  and  momentum  spaces  — the  variables  of  each  factor  in  the  split-operator 
propagator  were  chosen  such  that  the  Hamiltonian  in  each  factor  was  diagonal,  and  the 
coordinate  and  momentum  spaces  were  related  by  a Fourier  transformation.  A fast 
Fourier  transform  algorithm  was  employed  to  construct  the  kinetic  energy  propagator. 
A computer  code  in  Fortran  77  has  been  written  to  efficiently  calculate  femtosecond 
photodesorption  dynamics.  The  flow  diagram  for  the  wavepacket  propagation  computer 
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code  is  presented  in  Appendix.  This  program  package  has  provided  numerical  results 
for  CO  femtosecond  photodesorption  from  the  Cu(OOl)  surface. 

Results  and  Discussion 

The  calculations  of  the  CO  photodesorption  were  done  varying  several  different 
parameters:  the  light  pulse  fluence,  width,  and  central  or  carrying  frequency,  the 
magnitude  of  the  dissipative  imaginary  potentials  and  their  coupling,  and  the  temperature 
of  the  initial  system  and  the  isotope  masses  of  CO. 

Our  calculations  propagating  wavepackets  indicate  that,  for  indirect  desorption  by 
a light  pulse,  the  substrate  is  electronically  excited  and  deexcited  multiple  times  during 
the  desorption  and  leads  to  energy  transfer  to  nuclear  motions,  and  that  the  desorption 
dynamics  is  completed  in  times  comparable  to  the  pulse  width.  The  results  show  that  the 
desorption  yields  are  nonlinearly  dependent  on  light  pulse  fluence  which  is  consistent 
with  experimental  results.  The  desorption  happens  quicker  with  higher  frequency  light 
pulses  which  make  sense  because  each  photon  then  has  more  energy.  The  desorption 
yield  is  higher  for  the  light  pulse  with  the  same  fluence  but  wider  pulse  width  which  can 
be  explained  by  the  fact  that  there  is  more  time  for  energy  transfer  between  electrons 
and  nuclei.  The  desorption  yields  decrease  with  increasing  magnitudes  of  dissipation 
potentials  and  their  coupling;  this  is  because  even  the  coupling  can  increase  the  energy 
transfer  between  the  primary  and  secondary  regions,  and  the  substrate  will  dissipate 
more  energy  from  the  primary  region. 

The  calculated  photodesorption  yields  show  a nonlinear  dependence  on  light  pulse 
fluence.  With  the  light  pulse  parameters  chosen  the  same  as  that  used  in  the  experiment 
by  Prybyla  et  al.  [29],  and  the  effective  transition  dipole  scaled  to  produce  experimental 
yields  at  two  values  of  light  pulse  fluence,  and  the  dissipation  times  chosen  to  be 
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Tg  = 5p.s  and  = 25/.s  from  experiment,  our  results  are  in  consistent  with  the  remaining 
experimental  yield  dependences  on  light  pulse  fluences. 

The  calculated  photodesorption  delay  time  show  good  agreement  with  the  exper- 
iment by  Prybyla  et  al.  [29],  which  give  a desorption  time  of  less  than  325  fs  for  a 
100  fs  laser  pulse.  In  our  calculation,  the  desorption  time  is  around  300  fs  for  the 
same  light  pulse. 

Comparing  the  results  when  using  the  primary  transition  dipole  and  the  effective 
primary  transition  dipole,  the  photodesorption  yields  increased  dramatically  in  the  later 
case  especially  when  the  light  pulse  was  strong.  This  means  that  the  substrate  charge 
density  fluctuations  during  the  incident  of  the  light  pulse  have  a strong  effect  on  primary 
region  charge  density  redistribution  and  rapidly  increase  the  probability  of  electronic 
transitions  to  the  excited  state.  This  suggests  that  experimental  measurements  of  yields 
could  provide  an  estimate  of  the  effective  dipole  in  the  primary  region,  and  a check  on 
our  theoretical  model  of  dipole  enhancement  due  to  substrate  electronic  excitations. 

The  chirping  effects  on  photodesorption  were  also  calculated  and  shows  that  negative 
and  positive  chirping  have  quite  different  effects  on  CO  photodesorption  yields;  the 
negative  chirping  increase  the  CO  photodesorption  yields  and  the  positive  chirping 
decrease  CO  photodesorption  yields,  with  the  amount  of  changes  (increasing  and 
decreasing)  nonlinearly  dependent  on  chirping  strength.  The  chirping  alsosimilarly 
affects  the  population  of  electronic  ground  and  excited  states,  with  negative  chirping 
increasing  the  electronic  excited  state  population  and  positive  chirping  decreasing  it 

The  theory  developed  for  femtosecond  photodesorption  introduced  the  effective 
transition  dipole  and  incorporated  the  energy  dissipation  in  the  substrate  into  the  quantum 
dynamics  of  the  desorbing  species.  The  next  step  should  be  a detailed  study  of  the 
adsorbate-substrate  coupling  to  derive  expressions  for  the  effective  transition  dipole  and 
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dissipative  potentials  and  couplings  to  substitute  their  semiempirical  parametrization 
in  this  work.  Also,  the  theory  would  benefit  from  higher  quality  calculations  of  the 
excited  electronic  state  PES  to  improve  its  parametrization  and  reproduce  more  closely 
the  experimental  results. 

A Fortran  77  code  for  the  wavepacket  propagation  in  two  degrees  of  freedom  over 
two  coupled  complex  PESs  was  developed  and  applied  to  the  case  of  CO  undergoing 
frustrated — translation  vibrational  motions  and  desorbing  from  the  Cu(OOl)  surface. 
It  allows  for  a quantum  mechanical  calculation  of  any  two  dimensional  desorption 
process  described  with  only  Cartesian  coordinates.  However,  the  rotational  motion  of 
the  desorbing  species  is  not  taken  into  account  in  the  code  although  we  have  developed 
the  numerical  algorithm  by  combining  the  split-operator  scheme  with  a discrete  variable 
representation.  A further  improvement  of  the  code  would  be  to  implement  this  part  in 
the  code.  Also,  for  an  very  strong  light  pulse,  we  may  have  to  include  more  than  one 
excited  electronic  state  of  the  adsorbate-substrate  complex  involved  in  the  dynamics, 
so  the  algorithm  should  be  modified  to  include  another  PES  coupled  diabatically  to 
the  existing  ones.  This  would  obviously  increase  the  size  of  matrices  involved  in  the 
computation  and  also  make  the  transformation  between  the  electronic  representations 
more  complicated.  It  would  allow  description  of  additional  systems  although  the 
treatment  would  be  conceptually  unchanged. 


FLOW  CHART  OF  THE  PROGRAM  FOR  PROPAGATION  WITH  A 
TIME  DEPENDENT  HAMILTONIAN  OVER  TWO  COMPLEX 
POTENTIAL  ENERGY  SURFACES  IN  TWO  DEGREES  OF  FREEDOM 


The  program  (“photod.f’)  consists  of  four  functional  parts:  constructing  the  initial 
wavepacket  ( subroutine  “iwavepacket.f’),  constructing  potential  energy  surfaces  and 
their  couplings  (“pes.f ’),  propagating  the  wavepacket  one  time  step  (“propagate.f ’)  and 
calculating  properties  (“property .f’). 

The  subroutine  “iwavepacket.f’  constructing  the  initial  wavefunction  calls  the  fol- 
lowing subroutines: 

propdata.f  — read  in  the  propagation  data  file. 

propdata.IN  — input  data  for  propagation. 

The  subroutine  “pes.f’  constructing  the  potential  energy  surfaces  calls  the  following 
subroutines: 

pesfldata.f  — read  in  the  CONi(OOl)  potential  energy  surfaces  data  file. 

pesfldata.IN  — input  data  for  CONi(OOl)  photodesorption. 

pesf2data.f  — read  in  the  COCu(OOl)  potential  energy  surfaces  data  file. 

pesf2data.IN  — input  data  for  COCu(OOl)  photodesorption. 

The  subroutine  “propagate.f’  propagating  the  wavefunction  one  time  step  calls  the 
following  subroutines: 

adcp2d.f  — advances  a wavepacket  in  the  adiabatic  representation  by  half  a time 
step,  with  the  imaginary  potential  part  of  the  time  evolution  operator. 

adft2d.f  — advances  a wavepacket  one  time  step  with  the  kinetic  energy  operator 
in  momentum  space  using  a fast  Fourier  transform  (“four2d.f ’). 
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adpo2d.f  — advances  a wavepacket  by  half  a time  step,  with  the  real  potential  part 
of  the  time  evolution  operator. 

repr2d.f  — executes  a matrix  transformation  for  a basis  of  two  electronic  states. 
repi2d.f  — executes  an  inverse  transformation  for  a basis  of  two  electronic  states. 
four2d.f  — fast  Fourier  transform  between  coordinate  space  and  momentum  space. 
The  subroutine  “property .f’  calculating  the  time  dependent  physical  properties  calls 
the  following  subroutines: 

popu2d.f  — computes  the  populations  of  the  two  electronic  states  by  integrating  the 
wavepacket  amplitudes  at  given  time. 

yield2d.f  — computes  the  photodesorption  yields. 

acrf2d.f  — computes  the  autocorrelation  function  of  a wavepacket. 

The  calculation  results  were  included  in  the  following  data  files. 

YIELD. OUT  — desorption  yields  verse  time. 

POP  LOUT  — ground  population  distribution  verses  time. 

POP2. OUT  — excited  population  distribution  verses  time. 

AMPGR#.OUT  — ground  total  amplitude  distribution  in  certain  time. 
AMPER#.OUT  — excited  total  amplitude  distribution  in  certain  time. 


129 


f START 

^ Program  photod.f 


1 

Time-dependent 
Properties  (property.!) 


BIBLIOGRAPHY 


1.  Desorption  Induced  by  Electronic  Transitions  DIET  /,  Springer- Verlag,  New  York, 
1983,  N.  H.  Tolk,  M.  M.  Traum,  J.  C.  Tully  and  T.  E.  Madey,  editor. 

2.  M.  C.  Lin  and  G.  Ertl,  Ann.  Rev.  Phys.  Chem.  37,  587  (1986). 

3.  P.  Avouris  and  R.  E.  Walkup,  Ann.  Rev.  Phys.  Chem.  40,  173  (1989). 

4.  D.  S.  King  and  R.  R.  Cavanagh,  Adv.  Chem.  Phys.  76,  45  (1989). 

5.  R.  R.  Cavanagh,  D.  S.  King,  J.  C.  Stephenson,  and  T.  F.  Heinz,  J.  Phys.  Chem.  97, 
786  (1993). 

6.  F.  M.  Zimmerman  and  W.  Ho,  Surf.  Sci.  Reports  22,  128  (1995). 

7.  J.  Manz,  Molecular  wavepacket  dynamics:  Theory  for  experiments  1926-1996,  World 
Scientific,  Singapore,  in  the  press,  V.  Sundstrom,  editor. 

8.  D.  Menzel  and  R.  Gomer,  J.  Chem.  Phys.  40,  1164  (1964). 

9.  D.  Menzel  and  R.  Gomer,  J.  Chem.  Phys.  41,  3311  (1964). 

10.  P.  A.  Redhead,  Can.  J.  Phys.  42,  886  (1964). 

11.  P.  R.  Antoniewicz,  Phys.  Rev.  B 21,  3811  (1980). 

12.  J.  W.  Gadzuk,  L.  J.  Richter,  S.  A.  Buntin,  D.  S.  King,  and  R.  R.  Cavanagh,  Surf. 
Sci.  235,  317  (1990). 

13.  J.  W.  Gadzuk,  Phys.  Rev.  B 44,  13466  (1991). 

14.  J.  W.  Gadzuk,  Femtosecond  Chemistry,  pages  vol.2,  p.603,  VCH,  New  York,  1995, 
J.  Manz  and  L.  Woeste,  editor. 

15.  D.  M.  Newns,  T.  F.  Heinz,  and  J.  A.  Misewich,  Progr.  Theor.  Phys.  Suppl.  106, 
411  (1991). 


130 


131 


16.  J.  A.  Misewich,  T.  F.  Heinz,  and  D.  Newns,  Phys.  Rev.  Lett.  68,  3737  (1992). 

17.  M.  Brandbyge,  P.  Hedegard,  T.  F.  Heinz,  J.  A.  Misewich,  and  D.  M.  Newns,  Phys. 
Rev.  B 52,  6042  (1995). 

18.  B.  Fain  and  S.  H.  Lin,  Chem.  Phys.  161,  515  (1992). 

19.  B.  Fain,  S.  H.  Lin,  J.  Grotemeyer,  and  E.  W.  Schlag,  Chem.  Phys.  Lett.  302,  357 
(1993). 

20.  J.  C.  Tully,  M.  Gomez,  and  M.  Head-Gordon,  J.  Vac.  Sci.  Technol.  A , 1914  (1993). 

21.  C.  Springer,  M.  Head-Gordon,  and  J.  C.  Tully,  Surf.  Sci.  320,  L57  (1994). 

22.  M.  Head-Gordon  and  J.  C.  Tully,  J.  Chem.  Phys.  103,  10137  (1995). 

23.  F.  M.  Zimmerman  and  W.  Ho,  J.  Chem.  Phys.  100,  7700  (1994). 

24.  P.  Saalfrank,  R.  Baer,  and  R.  Kosloff,  Chem.  Phys.  Lett.  230,  463  (1994). 

25.  P.  Saalfrank  and  R.  Kosloff,  J.  Chem.  Phys.  105,  2441  (1996). 

26.  H.  Guo,  Chem.  Phys.  Lett.  240,  393  (1995). 

27.  L.  Liu,  H.  Guo,  and  T.  Seideman,  J.  Chem.  Phys.  104,  8757  (1996). 

28.  J.  A.  Prybyla,  T.  F.  Heinz,  J.  Misewich,  M.  M.  T.  Loy,  and  J.  H.  Glownia,  Phys. 
Rev.  Lett.  64,  1537  (1990). 

29.  J.  A.  Prybyla,  H.  W.  K.  Tom,  and  G.  Aumiller,  Phys.  Rev.  Lett.  68,  503  (1992). 

30.  F.  Budde,  T.  F.  Heinz,  M.  M.  T.  Loy,  J.  A.  Misewich,  and  F.  D.  Rougemont,  Phys. 
Rev.  Lett.  66,  3024  (1991). 

31.  F.  Budde,  T.  F.  Heinz,  A.  Kalamarides,  M.  M.  T.  Loy,  and  J.  A.  Misewich,  Surf. 
Sci.  283,  143  (1993). 

32.  J.  A.  Misewich,  A.  Kalamarides,  T.  F.  Heinz,  U.  Hofer,  and  M.  M.  T.  Loy,  J. 
Chem.  Phys.  100,  736  (1993). 

33.  T.  A.  Germer,  J.  C.  Stephenson,  E.  J.  Heilweil,  and  R.  R.  Cavanagh,  J.  Chem. 
Phys.  101,  1 (1994). 


132 


34.  F.  Kao,  D.  Busch,  D.  Cohen,  D.  G.  da  Costa,  and  W.  Ho,  Phys.  Rev.  Lett.  71, 
2094  (1993). 

35.  E.  Hasselbrink,  Laser  Spectroscopy  and  Photochemistry  on  Metal  Surfaces,  World 
Scientific,  Singapore,  in  the  press,  H.  L.  Dai  and  W.  Ho,  editor. 

36.  T.  Hertel,  M.  Wolf,  and  G.  Ertl,  J.  Chem.  Phys.  102,  3414  (1995). 

37.  D.  G.  Busch,  S.  Gao,  R.  A.  Pelak,  M.  F.  Booth,  and  W.  Ho,  Phys.  Rev.  Lett.  75, 

673  (1995). 

38.  D.  G.  Busch  and  W.  Ho,  Phys.  Rev.  Lett.  77,  1338  (1996). 

39.  L.  M.  Struck,  L.  J.  Richter,  S.  A.  Buntin,  R.  R.  Cavanagh,  and  J.  C.  Stephenson, 
Phys.  Rev.  Lett.  77,  4576  (1996). 

40.  D.  A.  Micha  and  D.  Srivastava,  Chem.  Phys.  Lett.  162,  376  (1989). 

41.  D.  Srivastava  and  D.  A.  Micha,  J.  Chem.  Phys.  94,  4900  (1991). 

42.  D.  Srivastava  and  D.  A.  Micha,  J.  Chem.  Phys.  95,  380  (1991). 

43.  D.  Beksic  and  D.  A.  Micha,  J.  Chem.  Phys.  103,  3795  (1995). 

44.  Z.  G.  Yi,  D.  Beksic,  and  D.  A.  Micha,  Intern.  J.  Quantum  Chem.  64,  71  (1997). 

45.  D.  A.  Micha  and  Z.  G.  Yi,  J.  Chem.  Soc.,  Faraday  Trans.  93,  969  (1997). 

46.  J.  A.  Pople  and  D.  L.  Beveridge,  Approximate  Molecular  Orbital  Theory,  McGraw- 
Hill,  New  York,  1970. 

47.  M.  C.  Zemer,  Review  of  Computational  Chemistry,  vol.  II,  pages  313-366,  VCH 
Publishing,  New  York,  1991,  K.  B.  Lipkowitz  and  D.  B.  Boyd,  editor. 

48.  M.  C.  Zemer,  Metal-ligand  Interactions,  pages  493-531,  Kluwer,  Dordrecht,  1992, 
N.  Russo  and  D.R.  Salahub,  editor. 

49.  M.  D.  Feit,  J.  A.  Fleck,  and  A.  Steiger,  J.  Comput.  Phys.  47,  412  (1982). 

50.  F.  M.  Zimmermann  and  W.  Ho,  Surf.  Sci.  Reports  22,  130  (1995). 

51.  P.  Chu,  S.  A.  Buntin,  L.  J.  Richter,  and  R.  R.  Cavanagh,  J.  Chem.  Phys.  101,  2929 
(1994). 


133 


52.  A.  Peremans,  K.  Fukutani,  K.  Mase,  and  Y.  Murata,  Phys.  Rev.  B 47,  4135  (1993). 

53.  A.  Peremans,  K.  Fukutani,  K.  Mase,  and  Y.  Murata,  Surf.  Sci.  283,  189  (1993). 

54.  C.  Springer  and  M.  Head-Gordon,  Chem.  Phys.  205,  73  (1996). 

55.  F.  T.  Smith,  Phys.  Rev.  179,  111  (1969). 

56.  M.  Baer,  Chem.  Phys.  15,  49  (1979). 

57.  J.  A.  Olson  and  D.  A.  Micha,  Intern.  J.  Quantum  Chem.  22,  971  (1982). 

58.  F.  M.  Fernandez  and  D.  A.  Micha,  J.  Chem.  Phys.  98,  8173  (1992). 

59.  C.  Kittel,  Quantum  Theory  of  Solids,  John  Wiley  and  Sons,  New  York,  1963. 

60.  D.  A.  Micha  and  K.  Runge,  Phys.  Rev.  Lett.  50,  322  (1994). 

61.  K.  Blum,  Density  Matrix  Theory  and  Applications,  Plenum,  New  York,  1981. 

62.  R.  Kosloff,  J.  Phys.  Chem.  92,  2087  (1988). 

63.  C.  Leforestier,  R.  Bisseling,  C.  Cerjan,  M.  D.  Feit,  and  R.Friesner,  J.  Comput. 
Phys.  94,  59  (1991). 

64.  R.  Kosloff,  Ann.  Rev.  Phys.  Chem.  45,  145  (1994). 

65.  R.  Kosloff  and  H.  Tal-Ezer,  Chem.  Phys.  Lett.  127,  223  (1981). 

66.  H.  Tal-Ezer  and  R.  Kosloff,  J.  Chem.  Phys.  81,  3967  (1984). 

67.  Z.  H.  Zhang  and  D.  J.  Kouri,  Phys.  Rev.  A A34,  2687  (1986). 

68.  R.  C.  Mowrey  and  D.  J.  Kouri,  J.  Chem.  Phys.  84,  6466  (1986). 

69.  T.  N.  Truong,  J.  J.  Tanner,  P.  Bala,  J.  A.  McCammon,  and  D.  J.  Kouri,  J.  Chem. 
Phys.  96,  2077  (1991). 

70.  R.  Kosloff  and  D.  Kosloff,  J.  Chem.  Phys.  79,  1823  (1983). 

71.  R.  W.  Heather  and  H.  Metiu,  J.  Chem.  Phys.  91,  1596  (1989). 

72.  X.-P.  Jiang,  R.  Heather,  and  H.  Metiu,  J.  Chem.  Phys.  90,  2555  (1989). 


134 


73.  M.  Jacon,  O.  Atabek,  and  C.  LeForestier,  J.  Chem.  Phys.  91,  1585  (1989). 

74.  A.  D.  Bandrauk  and  H.  Shen,  Chem.  Phys.  Lett.  176,  428  (1991). 

75.  J.  A.  Fleck,  J.  R.  Morris,  and  M.  D.  Feit,  Appl.  Phys.  10,  129  (1976). 

76.  W.  H.  Press,  S.  A.  Teukolsky,  V.  T.  Venerling,  and  B.  P.  Flannery,  Numerical 
Recipes  in  Fortran,  Cambridge  Univ.  Press,  Cambridge,  1992. 

77.  J.  Lill,  G.  Parker,  and  J.  Light,  Chem.  Phys.  Lett.  89,  483  (1982). 

78.  R.  W.  Heather  and  J.  C.  Light,  J.  Chem.  Phys.  79,  147  (1983). 

79.  J.  Light,  I.  Hamilton,  and  J.  Lill,  J.  Chem.  Phys.  82,  1400  (1984). 

80.  J.  Sheng  and  J.  Z.  H.  Zhang,  J.  Chem.  Phys.  96,  3866  (1992). 

81.  J.  Dai,  J.  Sheng,  and  J.  Z.  H.  Zhang,  J.  Chem.  Phys.  101,  1555  (1994). 

82.  H.  Wei  and  J.  T.  Carrington,  J.  Chem.  Phys.  101,  1343  (1994). 

83.  D.  Beksic,  Electronically  diabatic  photodesorption  of  molecules  adsorbed  on  metal 
surfaces,  PhD  thesis.  University  of  Florida,  Department  of  Physics,  1994. 

84.  S.Andersson  and  J.  B.  Pendry,  Phys.  Rev.  Lett.  43,  363  (1979). 

85.  J.  C.  Tracy,  J.  Chem.  Phys.  56,  2748  (1972). 

86.  C.  J.  Hirschmugl,  G.  P.  Williams,  F.  M.  Hoffman,  and  Y.  J.  Chabal,  Phys.  Rev. 
Lett.  65,  480  (1990). 

87.  J.  Ellis,  G.  Witte,  and  J.  P.  Toennies,  J.  Chem.  Phys.  102,  5059  (1995). 

88.  J.  E.  Demuth  and  P.  Avouris,  Phys.  Rev.  Lett.  47,  61  (1981). 

89.  P.  Nordlander,  Phys.  Rev.  B 46,  2584  (1992). 

90.  R.  M.  Jaeger,  K.  Homann,  H.  Kuhlenbeck,  and  H.-J.  Freund,  Chem.  Phys.  Lett. 
203,  41  (1993). 


91.1.  Kinoshita  and  A.  Misu,  J.  Chem.  Phys.  102,  2970  (1995). 


135 


92.  R.  Gomer,  Desorption  induced  by  Electronic  Transitions  DIET  /,  Vol.  24  of  Springer 
Series  in  Chemical  Physics,  page  40,  Springer- Verlag,  New  York,  1983,  N.  H.  Tolk 
and  M.  M.  Traum  and  J.  C.  Tully  and  T.  E.  Madey,  editor. 

93.  T.  E.  Madey,  J.  J.  T.  Yates,  D.  A.  King,  and  C.  J.  Uhlaner,  J.  Chem.  Phys.  52, 
5215  (1970). 

94.  J.  E.  Houston  and  T.  E.  Madey,  Phys.  Rev.  B 26,  554  (1982). 

95.  A.  Szabo  and  J.  John  T.  Yates,  J.  Chem.  Phys.  102,  563  (1995). 

96.  H.  W.  K.  Tom,  Amer.  Phys.  Soc.  meeting,  Los  Angeles,  March  (1998). 


BIOGRAPHICAL  SKETCH 


Zhigang  Yi  was  bom  on  June  26'*’  1964  in  Wuhan,  China.  He  was  taken  care  of  by 
his  grandmother  in  Hunan  until  he  was  seven  years  old  due  to  the  Culture  Revolution 
in  China.  All  of  his  early  education  took  place  in  Wuhan.  In  the  Summer  of  1985  he 
received  his  B.S.  degree  in  electronics  and  information  engineering  science  from  Wuhan 
University.  The  same  year  he  began  his  graduate  education.  In  the  Summer  of  1988  he 
received  his  M.S.  degree  in  space  physics  from  the  same  university.  After  graduation, 
he  worked  as  a software  engineer  and  assistant  manager  at  the  Wuhan  Instrument  Co. 
which  was  the  fourth  largest  instmment  company  in  China. 

In  the  Fall  1993  he  came  to  the  USA,  and  enrolled  in  the  Department  of  Physics  at 
the  University  of  Florida  for  further  studies.  He  joined  the  Quantum  Theory  Project  in 
the  Summer  of  1994.  After  a successful  completion  of  the  Ph.D  compressive  qualifying 
examination  in  the  fall  of  1994  he  was  admitted  to  the  Ph.D  program  and  began  his 
research  on  computational  and  surface  physics. 

After  completing  his  Ph.D  studies  at  UF,  Zhigang  is  prepared  to  find  an  industry 
position  in  the  computer  simulation,  visualization  and  animation  field  and  to  enjoy  many 
future  years  with  his  wife  Yipeng  and  his  son  Hongda  (lulu). 


136 


I certify  that  I have  read  this  study  and  that  in  my  opinion  it  conforms  to  acceptable 
standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality,  as  a 
dissertation  for  the  degree  of  Doctor  of  Philosophy. 


David  A.  Micha, 
Professor  of  Physics 


I certify  that  I have  read  this  study  and  that  in  my  opinion  it  conforms  to  acceptable 
standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality,  as  a 
dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Philip  f. B rue  at 
Associate  Professor  of  Chemistry 


I certify  that  I have  read  this  study  and  that  in  my  opinion  it  conforms  to  acceptable 
standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality,  as  a 
dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Charles  F.  Hooper,  Jr. 
Professor  of  Physics 


I certify  that  I have  read  this  study  and  that  in  my  opinion  it  conforms  to  acceptable 
standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality,  as  a 
dissertation  for  the  degree  of  Doctor  of  Philosophy. 

Pradeep  Kumar 
Professor  of  Physics 


I certify  that  I have  read  this  study  and  that  in  my  opinion  it  conforms  to  acceptable 
standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality,  as  a 
dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Michael  C.  Zerner  Z/ 
Distinguished  Professor  of  Physics 


This  dissertation  was  submitted  to  the  Graduate  Faculty  of  the  Department  of 
Chemistry  in  the  College  of  Liberal  Arts  and  Sciences  and  to  the  Graduate  School 
and  was  accepted  as  partial  fulfillment  of  the  requirements  for  the  degree  of  Doctor 
of  Philosophy. 


December  1998 

Dean,  Graduate  School 


